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ABSTRACT 


The  concept  of  transporting  solids  in  the  form  of  capsules 
in  pipelines  evolved  at  the  Research  Council  of  Alberta  during  the 
fifties.  Since  then  extensive  experimental  data  have  been  generated 
but  as  yet  no  satisfactory  correlations  are  available  covering  all 
the  physical  variables. 

The  data  available  from  the  Research  Council  of  Alberta  have 
been  analysed  using  Bayesian  statistical  methods.  Capacity 
correlations  have  been  obtained  for  pipe  sizes  ranging  from  1/2" 
to  10"  in  diameter.  Little  success,  however,  has  been  obtained  in 
modeling  the  pressure  drop  in  the  system. 

The  second  approach  to  the  problem  is  fluid  mechanics  modeling. 

A  preliminary  physical  model  envisaging  the  capsule  flow  as  consisting 
of  a  series  of  infinitesimal  parallel  slipper  bearings  is  given  here. 
Using  a  set  of  simplifying  assumptions  including  one  which  neglects 
secondary  flows,  capacity  and  pressure  drop  correlations  have  been 
developed  for  laminar  flow.  As  expected  from  the  assumptions  made, 
the  model  predictions  for  both  the  bulk  velocity  and  pressure  gradient 
for  a  given  capsule  velocity  are  lower  than  the  observed  values. 

However,  addition  of  only  one  empirical  parameter  in  each  correlation 
for  a  given  pipeline  system  gives  reasonable  agreement  with  the  observed 
values.  It  is  expected  that  substantial  improvements  in  predictions 
would  result  if  secondary  flows  can  be  incorporated  in  the  model. 
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1.  INTRODUCTION 


1.1  Capsule  Flow 

1.1.1  Historical 

The  development  of  solids  transportation  through  pipelines  has 
not  played  a  significant  role  in  our  overall  technological  growth.  The 
concept  was  suggested  in  California  as  early  as  1850.  The  first 
practical  application,  however,  did  not  materialize  until  1914  when  an 
8-inch  cast  iron  pipeline  1750  feet  long,  was  built  to  transport  coal 
slurries  from  the  London  docks  to  power  stations.  By  1959,  some  20 
pipelines  were  carrying  30  million  long  tons  annually  of  mixtures  of 
phosphate  rock,  sand  and  clay.  During  the  fifties  and  sixties,  pipelines 
were  installed  for  transporting  varieties  of  solids  including  ore 
concentrates,  limestone,  phosphate  and  coal  (1). 

New  developments  in  slurry  pipelining  depend  largely  on  solving 
a  host  of  technical  problems.  Some  of  the  problems  arise  in  processing 
the  solids  before  pumping,  keeping  the  solids  in  suspension  during 
transportation  and  separating  the  solids  from  the  carrying  fluid  and 
drying.  Abrasive  action  of  solids  in  the  pipeline  may  demand  special 
material  or  design  for  the  pipeline. 

A  novel  departure  from  the  established  solids  pipeline  methods 
was  suggested  by  the  findings  of  Charles,  Govier  and  Hodgson  (2)  on 
the  flow  of  oil  in  the  presence  of  water.  If  water  and  oil  of  the 
same  density  flow  through  a  pipeline,  stable  spherical  or  cylindrical 
slugs  of  oil  in  water  or  water  in  oil  can  be  formed.  Furthermore,  these 
bubbles  could  be  equal  density  solids  and  subsequent  experimental  tests 
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indeed  proved  that  this  mode  of  solids  transportation  was  feasible. 

The  idea  was  later  extended  to  solids  capsules  heavier  than  that  of 
the  carrying  liquid. 

1.1.2  Empirical  Work 

Experimental  work  has  been  going  on  in  the  area  of  capsule  flow 
at  the  Research  Council  of  Alberta  since  1958.  The  aim  of  the  program 
is  to  obtain  sufficient  understanding  of  the  capsule  flow  system  for 
a  rational  assessment  and  design  of  a  commercial  installation.  In  an 
attempt  to  obtain  the  dependence  of  capsule  velocity  and  pressure 
gradient  on  the  capsule,  the  carrying  liquid  and  the  pipe  characteristics, 
several  empirical  studies  have  been  reported  (3,  4,  5,  6,  7,  8,  9). 

The  bulk  of  the  work  has  been  done  at  the  Research  Council  of  Alberta, 
Edmonton,  Alberta.  Some  quantitative  relationships  have  been  obtained 
for  specific  conditions.  The  main  results  for  cylindrical  capsules 
may  be  summarized  as  follows. 

The  capsule  travels  at  a  slight  angle  to  the  axis  of  the  pipe, 
generally  with  a  nose-up  position.  The  clearance  at  the  capsule  bottom 
is  extremely  small,  being  barely  visible,  even  for  capsules  of  density 
nearly  equal  to  that  of  the  carrying  liquid  (say  within  2%) .  Over 
a  rather  wide  range  of  operating  conditions,  capsule  velocity  is  found 
to  be  roughly  equal  to  (within  +  20%)  the  bulk  velocity  of  the  liquid. 

The  capsule  length/diameter  ratio  and  end  shapes  have  little  effect  on 
the  capsule  velocity  and  pressure  gradients  for  capsule  length/diameter 
ratios  greater  than  6. 

The  capsule  velocity  is  also  found  to  vary  linearly  with  liquid 
bulk  velocity  except  at  very  low  velocities,  for  a  given  system  of 
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capsule,  carrying  liquid  and  pipeline.  Pressure  gradients,  however, 
follow  no  simple  relationship  with  the  liquid  bulk  velocity.  Velocity 
measurements  have  relatively  small  variance  and  are  generally  reproducible. 
While  pressure  drop  measurements  are  reproducible  on  the  run  to  run  basis 
it  was  found  they  were  not  reproducible  on  certain  pipelines  after  exten¬ 
sive  experimentation  suggesting  some  change  in  the  experimental  condi¬ 
tions.  Variances  of  the  capsule  velocity  measurements  and  that  of  the 
pressure  drop  measurements  vary  inversely  with  the  viscosity  of  the 
carrying  liquid,  with  the  capsule  to  pipe  diameter  ratio  and  with  the 
difference  between  the  capsule  and  liquid  densities.  These  inferences 
are  drawn  from  the  plots  of  the  raw  data  obtained  from  the  Research 
Council  of  Alberta. 

The  pressure  drop  is  lowest  and  the  velocity  ratio  (defined  as 
capsule  velocity  divided  by  liquid  bulk  velocity)  is  highest  when  the 
capsule  density  is  equal  to  that  of  the  liquid.  In  particular,  pressure 
drop  is  very  sensitive  to  capsule  density,  increasing  rapidly  with 
increase  in  capsule  density.  No  simple  relationship  can  be  established 
between  the  pressure  drop  and  capsule/pipe  diameter  ratio,  although, 
generally  the  velocity  ratio  increases  with  increasing  diameter  ratio. 

The  data  are  not  sufficient  to  establish  the  effect  of  the 
liquid  viscosity  on  the  pressure  drops  and  the  velocity  ratios. 

However,  there  are  indications  that  for  a  given  capsule/pipe  system, 
increasing  the  liquid  viscosity  from  a  low  initial  value,  decreases 
the  pressure  drop  while  velocity  ratio  increases  up  to  a  certain  point. 
Further  increases  in  viscosity  cause  an  increase  in  the  pressure  drop 
and  a  decrease  in  the  velocity  ratios. 
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Velocity  ratios  and  pressure  gradients  are  not  significantly 
different  for  a  train  of  capsules  carrying  many  capsules  from  those  for 
a  single  capsule.  The  line  fill  and  spacing  among  the  capsules  also 
has  only  minor  effects  on  the  velocity  ratios  and  pressure  drops. 

Another  interesting  empirical  observation  recently  made  by 
Ellis  and  Kruyer  (10)  is  that  pressure  gradients  are  substantially 
(75-80%)  reduced,  with  attendant  increase  of  the  capsule/liquid 
velocity  ratio,  by  fitting  the  front  half  of  the  capsule  with  a  collar. 
Sledge-type  runners  to  lift  the  capsule  off  the  bottom  of  the  pipe 
are  also  successful  in  reducing  the  pressure  drop  but  are  not  quite  as 
effective  as  the  collars. 

Some  attention  has  been  given  to  the  type  of  materials  that  can 
be  transported  and  the  techniques  for  transforming  solids  to  the  capsule 
form.  Materials  such  as  wheat,  chemicals  and  others  that  must  not  come 
in  contact  with  the  carrying  liquid  can  be  encapsulated  in  thin  metal  or 
plastic  containers.  Minerals  can  be  cast  in  the  form  of  solid  capsules. 
Yet  another  mode  is  to  slightly  extrude  solid  paste  to  form  cylindrical 
capsules  and  transport  these  so-called  paste  slugs  in  a  carrying  fluid 
that  is  of  substantially  different  nature  than  the  liquid  used  for 
making  the  paste.  An  example  is  a  paste  slug  formed  from  pulverized 
coal  and  oil  and  transported  with  water  as  the  carrying  medium.  The 
high  interfacial  tension  between  water  and  oil  keeps  the  slug  together 
and  this  mode  of  transportation  has  also  been  found  to  have  several 
other  operational  advantages  (11) . 
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1.1.3  Theoretical  Developments 

The  experimental  work  at  the  Research  Council  of  Alberta  has 
stimulated  some  theoretical  analysis  of  the  capsule  flow.  Without 
exception,  all  the  work  has  been  focussed  on  the  cylindrical  capsule 
with  an  effective  infinite  length.  Charles  (12)  and  Kennedy  (13) 
developed  theoretical  relationships  for  the  velocity  ratio  and 
pressure  gradient  for  a  concentric  flow  of  cylindrical  capsules.  Kruyer , 
Redberger  and  Ellis  (14)  carried  out  similar  analysis  for  laminar  ec¬ 
centric  flow  of  capsule  for  a  known  clearance  between  the  capsule  bottom 
and  the  pipe  wall. 

Garg  and  Round  (15)  suggested  that  the  bottom  clearance  could  be 
obtained  from  a  force  balance  on  the  capsule  if  the  pressure  gradient 
is  known.  The  procedure  rests  upon  the  supposition  that  the  flow 
past  a  capsule  can  be  described  as  the  sum  of  pressure  flow  through  the 
annulus  due  to  the  pressure  gradient  and  Couette  flow  due  to  the  drag 
force  of  the  capsule  wall.  The  method  is  applicable  for  both  laminar 
and  turbulent  flow  but  its  utility  is  limited  since  it  involves  complex 
numerical  integration  at  least  for  the  laminar  flow.  Moreover,  they 
assume  that  the  eccentricity  is  known.  The  reverse  method  for  predicting 
pressure  gradients  from  known  veiocity  ratios  has  not  been  found  to  be 
very  successful.  The  approach  has  been  applied  to  the  data  of  the 
Research  Council  of  Alberta  by  Kruyer  and  Garg  (16)  with  limited 
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The  problem  of  capsule  flow  is  exceedingly  complex  and  does  not 
seem  to  be  amenable  to  a  simple  solution.  However,  an  attempt  to 
solve  the  problem  in  an  integrated  fashion  utilizing  all  that  we  know 
so  far  about  the  capsule  flow  is  expected  to  lead  to  some  useful 
insights.  An  attempt  on  these  lines  is  made  in  Chapter  3. 

1.2  Experimental  Set-Up 

The  experimental  set-up  used  by  the  Research  Council  of  Alberta 
is  described  fully  in  references  (16)  and  (17).  The  exact  arrangement 
varied  somewhat  with  the  different  pipelines  used  to  generate  the 
data.  Figure  1.1  shows  a  typical  schematic  representation  of  a 
horizontal  straight  shuttle  pipeline  that  allows  the  shuttling  back  and 
forth,  of  cylinderical  capsules.  Measurements  are  performed  when  the 
capsule  travels  through  the  test  section.  The  capsule  is  inserted 
in  a  lock  between  valve  2  and  3  and  is  held  stationary  until  the  liquid 
attains  steady  state.  The  capsule  is  then  injected  into  the  flowing 
liquid  by  closing  valve  3  and  opening  valve  2.  It  gathers  momentum 
rapidly  and  only  a  short  distance  is  needed  to  attain  steady  state 
capsule  velocity.  A  four-way  valve  (V^)  determines  the  direction  of 
liquid  flow  and  two  throttle  valves  (V^.,  V^)  are  used  to  control  the 
flow  rate. 

Up  to  three  centrifugal  pumps  can  be  used  in  parallel  to  provide 
a  flow  of  100  USGPM  or  liquid  bulk  velocities  up  to  8  feet/second. 
Pressure  sensitive  switches  are  provided  for  safety,  generally  set 
to  turn  off  the  pumps  if  the  pressure  exceeded  60  psi.  Sensitive 
photo-cells  also  turned  off  the  pumps  before  the  capsule  reached  the 
pipe  end  to  avoid  shock  to  the  pumps  of  sudden  capsule  stopping.  A 
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cooling  unit  connected  to  a  heat  exchanger  in  the  injection  tank 
maintains  temperature  control  within  2°F. 

Variations  in  the  internal  pipe  diameter  in  the  test  section 
can  be  kept  within  0.01  in.  except  for  the  10  in.  line  in  which  the 
maximum  variation  tolerated  is  0.03  in.  Capsule  length  can  be  kept 
constant  within  0.05  in.  and  maximum  variation  in  the  capsule  diameter 
can  be  held  at  0.03  in.  Capsule  volumes  and  weight  measurements  allow 
the  capsule  densities  to  approach  the  selected  densities  within  0.5%. 

A  PDP-8-S  (Digital  Equipment  Corporation)  mini-computer  is  used 
for  acquisition  and  partial  processing  of  data.  Figure  1.2  shows  the 
flow  of  data  and  control  information  between  the  central  processor  and 
the  equipment.  The  numbers  near  the  top  of  the  figure  indicate  how 
many  channels  were  used  to  record  the  variable. 

Temperatures  are  measured  at  the  pipe  inlet  and  outlet  by 
thermistors  that  are  regulated  by  a  1500  volt  supply.  The  digitized 
voltage  signal  from  the  thermistors  are  fed  to  binary  counters  that 
are  interfaced  with  the  central  processor.  Pressure  transducers  used 
to  measure  pressure  differentials  also  operate  on  a  similar  mechanism. 
The  voltage  from  one  of  the  pressure  transducers  can  also  be  displayed 
on  an  oscilloscope  to  allow  monitoring  of  the  pressure  differential. 

The  liquid  flow  is  monitored  by  a  magnetic  flow  meter  and  photo¬ 
cells  are  employed  to  detect  the  capsule.  A  crystal  oscillator  and  a 
counter,  interfaced  with  the  central  processor,  are  used  to  record  time. 
Capsule  velocities  are  measured  to  within  0.3%  over  the  entire  range 
and  liquid  flow  rates  are  within  1%.  The  pressure  transducers  allow 
measurements  to  within  0.2%  of  the  transducer  range  during  calibration. 
Over  a  period  of  a  month,  however,  the  calibration  error  can  be  as 


much  as  0.5%. 
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A  typical  trace  of  pressure  drop  between  the  pipe  taps  as  a 
capsule  moves  through  the  pipeline  is  shown  in  Figure  1.3.  In  a  typical 
pass  the  capsule  initially  moves  in  the  pipeline  towards  the  first 
pressure  tap  so  that  the  pressure  differential  is  caused  by  the 
flowing  liquid  only.  The  trace  of  the  pressure  differential  rises  as 
the  capsule  moves  past  the  first  tap  and  becomes  constant  as  the 
entire  capsule  length  comes  within  the  two  pressure  taps.  The  decline 
in  the  trace  corresponds  to  the  exit  of  capsule  past  the  second  tap. 

The  average  pressure  differential  is  computed  by  numerical  integration 
from  the  point  D  to  H  on  the  trace,  i.e. 

H 

AP  dt 

J 

Average  AP  =  AP  =  tj -  (1*1) 

m  n 

dt 

y 

D 

where  AP  is  the  instantaneous  pressure  differential  between  the  pressure 
taps  and  t  is  time.  A  similar  integration  is  carried  out  for  all  the 
other  variables.  The  pressure  gradient  for  the  capsule  PGC  can  be 
obtained  from  the  measured  average  pressure  differential  AP^  by 
the  following  formula: 

AP 

PGC  =  [AP  -  (L  -  L  )  -j~ ]/L  (1*2) 

m  P  c  J-'p  u 

where  L  is  the  pipe  length  between  the  pressure  taps,  L  is  the  capsule 
P 

length  and  AP^  is  the  pressure  differential  for  liquid  flow  alone. 
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1. 3  Model  Discrimination  and  Model  Building 
1.3.1  Mathematical  Modeling 

A  mathematical  model  is  an  abstract  formulation  aimed  to  simulate 
a  real  world  phenomenon.  A  true  mathematical  model  accounting  completely 
for  the  physical  phenomenon  is  a  hypothetical  concept  and  is  seldom 
realized  in  non-trivial  situations.  Nevertheless,  it  is  a  useful 
tool  for  the  scientists  and  engineers  in  representing  the  important 
characteristics  of  the  actual  methanism.  One  is  often  concerned  with 
a  situation  where  a  dependent  variable  y,  is  related  to  a  set  of 
independent  variables,  x  =  (x^,  x^ ,  ...  x^}.  For  example,  an  engineer 
might  be  studying  the  dependency  of  the  pressure  gradient  y  through  a 
pipe  on  the  pipe  diameter  x^,  the  liquid  viscosity  x the  liquid 
density  x^,  and  others.  Theoretically  there  exists  a  relation  of  the 
type 

y  =  <J>(x,  jO  (1*3) 

between  the  dependent  variable  y  and  the  set  of  independent  variables 
x,  where  a  =  {a±i  a2,  ...  ar>  is  the  set  of  r  parameters  which  may  be 
the  physical  constants  of  the  system  under  study.  Usually,  a  simpler 
relation  that  closely  approximates  the  outstanding  features  of  (1.3) 
is  satisfactory. 

In  the  course  of  the  modeling  process,  the  analyst  is  faced  with 
two  classes  of  problems.  One  is  an  engineering  situation  where  the 
purpose  is  to  obtain  a  quantitative  description  of  the  process  for 
design  when  extensive  data  are  available  in  the  region  of  interest. 

The  "response  surface"  methodology  deals  with  this  class  of  problems 
and  an  empirical  polynomial  fit  is  usually  satisfactory  (18,  19,  20). 
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The  other  type  of  situation  arises  when  the  purpose  is  to 
establish  the  true  nature  of  the  process  and  several  mechanistic  models 
can  be  postulated.  It  is  desired  to  discriminate  among  these  rival 
mechanisms  so  as  to  identify  the  true  mechanism.  Here,  the  problems 
are  twofold.  In  the  first  place.  Informative  experiments  have  to  be 
designed;  secondly,  the  resulting  data  have  to  be  analysed  to  discriminate 
among  the  rival  model  forms. 

'  In  the  methodology  to  follow  in  the  next  chapter,  a  blend  of 
the  above  two  approaches  has  been  employed.  For  the  capsule  flow  problem, 
enough  mechanistic  insight  is  not  available  to  warrant  the  later  approach. 
Again,  a  purely  empirical  approach  would  require  long  and  unwieldy 
polynomials  since  the  number  of  independent  variables  is  rather  large. 
Semi-empirical  models  are  generated  on  the  basis  of  physical  reasoning 
or  merely  intuitive  judgement.  In  addition,  some  statistical  tools 
are  available  for  model  improvements  on  an  empirical  basis. 

1.3.2  Model  Discrimination  and  Parameter  Estimation 

Once  the  models  have  been  generated  either  empirically  or 

mechanistically,  discrimination  among  the  rival  model  forms  is  carried 

out  along  Bayesian  lines.  The  procedure  follows  that  developed  by 

D.  Quon  and  the  author  (21,  22).  Suppose  we  have  m  plausible  model 

forms  relating  the  dependent  variable  y  to  the  set  of  independent 

variables  x,  each  model  form  ^  containing  a  vector  of  parameters  a^ 

of  dimension  r.: 

3 

y  =  <J>  (x,  cO  j  =  1»  •••»  m  (l*^) 

j  “  ~~3 

and  n  sets  of  experimental  data,  i.e.  {y^,  2L^:  »  •••» 


■i 

- 

■ 


■ 
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If  there  Is  no  prior  information  about  the  parameters  in  the 


region  of  interest  then  the  prior  parameter  distribution  is  approximated 
by  the  improper  uniform  distribution  which  written  in  normal  form  is 
given  by 


(1.5) 


where  is  any  finite  vector  corresponding  to  the  mean  value  for 
aj  and  Cj ,  the  prior  variance  covariance  matrix,  has  infinite  values 


along  its  diagonal. 

A 

The  least  squares  estimate  is  found  by  solving  the  following 
extremization  problem: 


The  algorithm  for  the  parameter  estimation  is  simple  involving  only 
a  matrix  inversion.  Iterative  convergent  schemes  have  to  be 
employed  for  models  that  are  non-linear  with  respect  to  the  parameters. 
Marquardt’s  (23)  algorithm  has  been  used  for  non-linear  parameter 
estimation  in  this  work. 


The  unbiased  variance  of  fit  is  simply 


r 


j 


(1.7) 


th 

and  the  design  matrix  for  the  j  model 


U.  =  (u, , :  i=l ,  ...  n} 

“1  “ji 


(1.8) 
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where 


3  4> 


8a 


jk 


k  =  1,  ...  rJ 


jk 

x 


jk 

5i 


(1.9) 


The  information  matrix  is  given  by 


The  posterior  parameter  distribution  is  now: 


(1.10) 


f"(V  = 


(1.11) 


For  most  purposes,  it  is  adequate  to  assume  that  all  the  m  models 
are  equally  likely  before  any  experimental  observations  have  been  made. 
That  is 


- 


(1.12) 


The  posterior  probability  of  model  form  <j>^  is  obtained  by 
using  Bayes'  theorem 


p,(V 


P1  (fj)  LCzIfj) 


l 


k=l 


L<xl  V 


(1.13) 
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where 


L(zUj) 


L(llf.  >  f"(a^)  da 

h 


n  L(y,  <J>,  ,  jL, )  f"(0  d  a.  (1.14) 

J  1  J  “3  “3 

nfN<yihj(5i»  Sj),  °j}  d  3^ 


The  above  integral  has  been  evaluated  in  references  (21,  22) 
The  posterior  probability  p"(<J>j)  is  used  as  a  measure  of  the 
validity  of  the  model  form.  It  is  a  measure  of  the  probability  that 
the  given  model  is  the  true  model  out  of  the  set  of  models  under 
consideration. 


1.3.3  Tests  of  Model  Adequacy 

The  posterior  model  probability  is  only  a  measure  of  the 
relative  success  of  the  given  model  in  describing  the  data.  When  all 
the  perceived  models  are  poor  ones,  even  the  most  successful  may  not 
be  adequate  for  a  given  purpose.  Simple  tests  of  model  adequacy 
have  been  developed  in  references  (21)  and  (24)  and  may  be  summarized 
as  follows: 

a)  The  residual  errors  should  be  uncorrelated  with  all  the 

independent  variables.  Existence  of  any  such  correlation  suggests 
systematic  modeling  errors  and  that  the  model  can  be  further 


improved. 


v 
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b)  The  estimated  error  variance,  as  computed  by  relation  (1.7) 
should  not  be  larger  than  the  experimental  error  variance, 
estimated  on  physical  grounds.  •  The  variance  of  fit  should  not 
change  drastically  if  various  subsets  of  the  data  are  employed  or 
when  additional  data  are  incorporated. 

c)  The  parameter  estimates  should  not  vary  widely  with  the 
incorporation  of  additional  data.  Relatively  stable  parameter 
estimates  over  the  entire  operating  space  ensure  that  the  model 
under  consideration  is  general  and  not  merely  a  'local'  approxi¬ 
mation  of  some  other  intrinsic  model  form. 

1.3.4  Model  Improvement  and  Model  Generation 

If  all  the  perceived  models  are  found  to  be  inadequate  against 
the  criteria  of  section  (1.3.3),  the  models  need  to  be  improved  by  the 
addition  and/or  modification  of  some  terms.  Indeed,  entirely  new 
models  might  be  necessary.  However,  there  are  no  established  procedures 
available  for  model  generation  or  improvement.  In  dealing  with  the 
problem  of  the  capsule  flow,  it  was  found  that  many  of  the  models 
had  to  be  generated  on  an  ad  hoc  basis  and  most  improvements  were 
peculiar  to  the  situations  encountered.  Some  general  guidelines 
in  this  regard  that  could  be  isolated  from  the  experience  with  the 
capsule  flow  problem  are  given  below. 

a)  Examine  the  correlation  of  the  residual  errors  with  each  of  the 
independent  variables.  If  any  regular  trend  in  the  errors  is 
observed  while  examining  against  a  particular  variable  x^ ,  then 
the  model  can  be  improved  by  modifying,  adding  or  deleting  the 
terms  involving  x^ . 


■  .1.  I  y.  9  V 
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b)  The  marginal  variances  of  the  parameter  estimates  (the  diagonal 
elements  of  the  variance  covariance  matrix)  are  useful  in  model 
simplification.  A  large  parameter  variance  implies  that  the 
corresponding  term  may  be  safely  dropped  without  significant 
deterioration  in  the  model  performance. 

c)  Graphical  display  of  the  dependent  variables  versus  each  of  the 
independent  variables  can  sometimes  suggest  the  approximate 
functional  dependency.  Common  mathematical  functions  for  several 
values  of  the  parameters  may  be  plotted  on  the  same  scale  as  the 
actual  data. 

d)  Simple  transformation  of  the  independent  and  dependent  variables  is 
a  powerful  tool  for  generating  useful  function  groups.  Among 

the  common  transformations  are  logarithmic,  exponential,  power 
and  polynomial.  Graphical  display  of  the  dependent  variable  with 
all  other  independent  variables  held  constant,  if  such  data  are 
available,  is  usually  sufficient  to  suggest  a  suitable  transformation. 

e)  Semi-theoretical  physical  reasoning  can  help  generate  some  useful 
model  forms.  The  concept  of  driving  force  and  resistance  can  be 
employed  in  examining  the  functional  dependency.  Some  models  may 
be  discarded  as  being  incongruous  with  the  physics  of  the  problem. 
Several  others  may  be  generated  by  modification  of  known  models 
for  the  limiting  cases  (e.g.  that  of  the  concentric  cylinders). 

The  capsule  flow  modeling  problem  has  been  attempted  here  through 
two  different  approaches.  The  first  utilizes  some  of  the  theoretical 
insights  based  mainly  on  the  experimental  data  and  is  described  in 
Chapter  2.  The  alternate  attempt  while  incorporating  all  the 


' 
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important  empirical  observations  is  primarily  a  theoretical  analysis 
of  the  capsule  flow.  The  latter  constitutes  the  subject  matter  of 
Chapter  3.  The  last  chapter  sums  up  the  results  of  this  work. 
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2.  EMPIRICAL  CORRELATIONS 


2 . 1  Scope 

In  any  engineering  design  involving  the  transportation  of 
capsules  in  a  pipeline,  it  is  necessary  to  develop  quantitative 
relationships  between  such  key  design  variables  as  capsule  velocity 
(which  is  directly  related  to  pipeline  capacity)  and  pressure  drop  and 
other  parameters  such  as  pipeline  geometry,  fluid  properties,  and  the 
physical  characteristics  of  the  capsules.  Unfortunately,  the  complexities 
of  the  flow  phenomenon  have  frustrated  any  attempts  up  to  no  ft  to 
develop  a  general  theoretical  mechanistic  model  that  would  describe 
and  explain  capsule  pipeline  flow.  Consequently,  we  are  reduced  to 
formulating  empirical  correlations  based  on  the  extensive  but  by  no 
means  exhaustive  data  gethered  over  the  past  number  of  years,  in 
particular,  by  the  Research  Council  of  Alberta.  Difficulties  are 
encountered  even  in  developing  empirical  correlations.  Thus  while 
reasonably  satisfactory  capacity  relationships  will  be  presented  in 
this  chapter,  no  general  satisfactory  pressure  drop  correlation 
has  been  developed. 

2 . 2  Physical  Variables 

Even  a  cursory  physical  analysis  of  capsule  pipeline  flow  would 
indicate  that  at  least  the  following  independent  variables  might 
contribute  to  a  complete  specification  of  the  flow  system. 


, 
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D  pipe  diameter,  cm 

d  capsule  diameter,  cm 

capsule  length,  cm 

Sc  capsule  end  shape,  dimensionless  factor 

o  capsule  density,  g/cc 

p  liquid  density,  g/cc 

y  liquid  viscosity,  centipoise 

bulk  liquid  flow  velocity  in  the  pipe,  ft/sec 

Ep  roughness  factor  of  the  pipe  wall,  dimensionless 

Ec  roughness  factor  of  the  capsule  surface,  dimensionless 

Given  the  above,  the  following  may  be  considered  to  be  dependent 
or  derived  variables: 

capsule  velocity,  ft /sec 

k  ratio  of  capsule  diameter  to  the  pipe  diameter  d/D,  dimensionless 

V  annulus  velocity  (VK-k^V  )/(l-k^),  ft/sec 

sn  D  c 

VR  velocity  ratio  V  /V,  ,  dimensionless 

c  b 

PGF  pressure  gradient  when  liquid  flows  at  the  rate  of  V,  ,  psi 

b 

PGC  pressure  gradient  for  the  capsule  system,  psi 
PR  pressure  ratio  PGC/PGF,  dimensionless 

e  minimum  clearance  between  the  capsule  surface  and  pipe  wall,  cm 

n  coefficient  of  friction  between  the  capsule  surface  and  the 

pipe  wall,  dimensionless 

kinematic  viscosity  of  the  liquid  y/p,  centistokes 


v 
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2 . 3  Experimental  Data 

Obviously  not  all  the  independent  variables  could  be  thoroughly 
investigated,  nor  could  all  the  dependent  variables  be  measured.  In 
their  experimental  program,  except  for  certain  preliminary  studies, 
the  Research  Council  of  Alberta  concentrated  on  those  ranges  of 
parameter  values  which  would  most  likely  be  encountered  in 
industrial  applications.  The  experimental  procedure  has  been  described 
in  section  1.2.  The  ranges  of  conditions  under  which  the  Research 
Council  of  Alberta  obtained  data  are  summarized  in  Table  2.1. 

Extensive  runs  were  carried  out  for  each  different  type  of 
pipe  and  pipe  size  but  the  carrying  liquid  was  either  oil  or  water, 
thus  limiting  the  range  of  liquid  properties  investigated.  Both  the 
diameter  ratio  and  density  of  the  capsules  were  varied  over  substantial 
ranges.  All  of  the  independent  variables  except  Ep,  E^  and  were 
measured.  The  roughness  factors  Ep  and  E^  had  to  be  inferred; 
they  may  have  changed  during  the  course  of  experimentation,  and  such 
changes  may  indeed  have  been  responsible  for  the  apparent  change 
in  some  of  the  pressure  gradients  to  be  noted  later . 

Of  the  dependent  variables,  the  annulus  velocity  V  ,  the  clearance  e 
and  the  coefficient  of  friction  ri  could  not  be  measured. 

2. A  Capacity  Correlations 

It  was  attempted  to  establish  relationships  of  the  type 

Vc  =  f(Vb,  D,  k,  a,  p,  v,  Lc/d,  Ep,  Ec,  Sc)  (2.1) 

between  the  dependent  variable  and  the  set  of  independent  variables. 
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For  any  given  run,  all  the  independent  variables  except 
Vfa  have  been  held  constant.  The  first  step  towards  identifying  the 
relation  (2.1)  was  to  model  all  the  individual  runs  for  each  pipe 
size.  Among  others,  the  following  models  were  considered  to  fit  the 
data  of  individual  runs: 


Model  1 


VR  = 


ClVb 


d  +  c2vbJ) 


(2.2) 


Model  2 


V  «  b  +  b_V  ' 
c  1  2  b 


(2.3) 


where  c^,  c^,  c^,  b^,  and  b^  are  functions  of  D,  k,  a/p,  v,  iWd, 

E  ,  E  and  S  .  Since  these  independent  variables  have  constant  values 
pc  c 

for  a  given  run,  the  parameters  c^,  c2,  c^,  b^,  b2>  and  b^  can  be 
uniquely  estimated  for  each  run.  Model  1  was  formulated  from  the 


shapes  of  VR  versus  V  plots  and  model  2  was  inferred  from  V.  versus 

c  b 

V  plots  for  individual  runs. 


c 

Fitting  models  1  and  2  to  the  data  of  several  sample  runs 
resulted  in  the  lower  variance  of  fit  for  model  2  in  all  cases. 
Furthermore,  model  1  was  found  to  be  inadequate  when  tested  against  the 
criteria  of  section  1.3: 

1)  Widely  varying  estimates  of  parameters  depending  upon  the  region 
of  the  subset  of  data  used;  e.g.,  low  or  high  region. 

2)  Estimated  variance  of  fit  was  found  to  be  dependent  upon  the 
number  of  data  points  used  in  a  given  subset  and  the  region  of 


data  in  the  given  subset. 
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Model  2,  on  the  other  hand  resulted  in  relatively  stable 
parameters  and  variance  estimates.  Furthermore,  it  could  be  inferred 
from  parameter  estimates  of  only  a  few  sample  runs  that  within  95% 
confidence  limits,  b^  =  1.  Table  2.2  gives  the  values  of  the  parameter 
b^  and  its  95%  confidence  limits  for  some  sample  runs.  A  newly 
acquired  APL  plot  subroutine  (Appendix  G)  was  used  to  detect  any 
curvature  in  the  remaining  runs  and  fairly  good  straight  lines  were 
obtained.  This  obviated  the  need  for  testing  the  remaining  runs  by 
an  extensive  non-linear  search  for  the  parameter  b^.  Figure  2.1  shows 
the  plot  of  versus  for  a  typical  run. 

Model  2  was  thus  simplified  to 


V  =  b  +  b0V, 
c  1  2  b 

(2.4) 

where 

b..  —  f-  (k,  D,  a/p,  v,  L  /d,  E  ,  E  ,  S  ) 
i  i  c  p  c  c 

(2.5) 

and 

b.  =  f_(k,  D,  a/p,  v,  L  /d,  E  ,  E  ,  S  ) 

2  2  c  p  c  c 

(2.6) 

The  correlation  (2.4)  described  accurately  99%  of  the  experi¬ 
mental  data  within  +  1%  of  the  average  value  for  most  of  the  runs  and 
within  +  5%  for  all  the  runs.  Table  2.3  shows  the  parameter  and  variance 
estimates  from  the  data  on  the  1  1/4"  copper  pipeline. 

The  next  step  was  to  identify  the  functional  dependence  of  b^ 
and  b£  on  the  remaining  independent  variables.  A  few  preliminary 
runs  indicated  that  the  capsule  end  shape  factor  Sc  and  the  capsule 
length  to  diameter  ratio  iWd  did  not  materially  affect  the  parameters 

b.  and  b0  for  L  /d  greater  than  about  6  (5).  Unfortunately, 
12c 

the  effect  of  pipe  diameter  could  not  be  isolated  since  each  pipe 
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TABLE  2.2 

ESTIMATES  OF  FOR  SAMPLE  RUNS 
IN  I  1/4"  COPPER  PIPE  LINE 


Run  No.  95%  Limits 


428 

1.003 

+  .011 

429 

.989 

.010 

.2,  413 

.987 

.004 

421 

.977 

.007 

452 

.995 

.188 

447 

.929 

.231 

446 

1.073 

.106 

441 

.993 

.003 

440 

.994 

.004 

408 

1.002 

.003 

416 

1.005 

.003 

' 


* 

*  • 


Vc  (ft/sec) 
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Figure  2.1  Vc  versus  Vb  for  run  &  424,  I  1/4" 
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TABLE  2.3 

PARAMETERS  AND  VARIANCE  ESTIMATES 
FOR  1  1/4"  COPPER  PIPE  LINE 


Run  No. 

0 

k 

bi 

b2 

2 

a* 

408 

1.031 

.94 

-.03105 

1.05960 

1.91xl0"5 

410 

1.031 

.8895 

-.06898 

1.11423 

6.88xl0~5 

412+413 

1.032 

.839 

-.19789 

1.17697 

8.62xl0'5 

414 

1.031 

.7916 

-.27800 

1.22815 

1.93x10 

416 

1.125 

.94 

-.04387 

1.06434 

3.91xl0~5 

418 

1.125 

.8895 

-.19648 

1.13260 

8.82xl0~5 

420+421 

1.125 

.839 

-.45320 

1.2077 

2.63x10 

422 

1.125 

.7916 

-.73822 

1.27162 

3.38xl0~4 

424 

1.25 

.94 

-.09543 

1.06725 

6.31xl0"5 

426 

1.25 

.8895 

-.33449 

1.1457 

1.80xl0"4 

-4 

428+429 

1.25 

.839 

-.69178 

1.2175 

2.44x10 

430 

1.25 

.7916 

-1.44900 

1.3066 

2.34xl0~3 

432 

1.5 

.94 

-.17152 

1.07468 

3.09xl0"5 

434 

1.5 

.8895 

-.51370 

1.1682 

6.41xl0~5 

436+437 

1.5 

.839 

-.89606 

1.23461 

4.94x10 

438 

1.5 

.7916 

-1*72840 

1.34691 

5.52x10 

5.56xl0"5 

A 

4401-441 

2.0 

.94 

-.30188 

1.08865 

442 

2 

.8895 

-.76826 

1.157783 

6.84x10 

444+445 

2 

.839 

-1.38076 

1.2117848 

7.16x10 

446+447 

2 

.7916 

-2.3261 

1.328178 

1.50x10 

448 

3 

.94 

-.41845 

1.08540 

1.15x10 

450 

3 

.8895 

-1.0493 

1.16786 

7.51x10  J 

452+453 

3 

.839 

-2.2445 

1.27304 

1.90x10 

454 

3 

.7916 

-3.51224 

1.44130 

3.79x10 

456 

5 

.94 

-.57372 

1.07729 

1. 44x10" 

458 

5 

.8895 

-1.85227 

1.18449 

1.20x10 

' 

> 
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diameter  was  associated  with  a  different  pipe  material  (with  variation 

in  E  )  or  with  a  different  capsule  material  (with  variation  in  E  ) . 

P  c 

Therefore,  the  best  that  could  be  done  was  to  obtain  the  functional 
forms  for  the  parameters  b^  and  b^  for  each  experimental  pipeline 
rather  than  a  general  correlation  covering  all  the  pipelines. 

Thus 

bl  =  ^p^k»  °^P*  (2.1) 

-  <j>1(k,  a/p)<f>2(v)  (2.8) 

and 

b?  =  (k,  a/p,  v)  (2.9) 

z.  p 

-  (k,  o/p)ip2(\>)  (2.10) 


The  functional  forms  4>^  and  now  characterize  the  specific 
pipeline  size,  the  pipe  material  and  the  capsule  material.  The 
factorization  of  the  functional  forms  in  (2.8)  and  (2.10)  implies 
independence  of  the  effects  due  to  capsule/pipe  diameter  ratio  and 
capsule  density  on  the  one  hand  and  liquid  viscosity  on  the  other. 
There  is  no  theoretical  justification  for  the  factorization  but  for 
the  two  pipelines  in  which  viscosity  was  varied,  the  use  of  the 
empirically  best  forms  for  <f>2  and  \p2. 


JL _ 

0.25 

v 


(2.11) 


2 


=  v 


(2.12a) 


or 


(2.12b) 


gave  satisfactory  correlations.  In  the  other  pipeline  data,  the  range 
of  viscosity  variation  was  not  sufficient  to  test  the  validity  of  the 


assumption. 


' 
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The  functions  <f>^  and  were  found  to  be  more  complicated 
and  no  single  model  was  successful  for  all  the  pipelines,  suggesting 
that  any  functional  forms  not  including  surface  roughness  could  only 
give  locally  correct  predictions.  The  technique  used  for  selecting 
models  and  estimating  parameters  follows  along  the  Bayesian  lines 
described  in  section  1.3.  For  preliminary  modeling  of  <J>  and 
the  parameter  estimates  b^  and  b£  from  individual  runs  were  used.  The 
model  formulation  and  evaluation  consisted  of  examination  of  plots, 
residual  errors,  variance  of  fit  and  application  of  logarithmic  and 
power  transformation  of  variables. 

To  keep  the  analysis  tractable  and  the  resulting  functional 
forms  simple,  only  single  term  functional  forms  were  investigated. 

One  of  the  powerful  tools  for  generating  model  forms  employed  was 
semi-theoretical  reasoning.  For  example  in  determining  the  functional 
form  of  b ^  for  laminar  flow: 


Resistance  to  flow  “ 


1 


(clearance  =  1-k) 


c 


Driving  force 


a  surface  area  of  capsule  k 

or  the  two  surface  areas  (1+k) 

2 

«  cross  sectional  area  k 


V, 


So  that 


V 


Driving  force 

Resistance 


(l-k)C  k2 

(1+k) 


V, 


- 


> 
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Thus 


b2  “ 


2  c 

d+k) 


where  c  is  an  arbitrary  power  to  be  determined  from  the  use  of  data. 

As  an  example.  Tables  2.4  and  2.5  give  some  of  the  model  forms 
tried  for  and  \p^  using  parameter  estimates  b^  and  b 2  from  the  data 
of  1  1/4"  copper  line  carrying  steel  capsules.  Some  of  the  more 
successful  forms  for  <f>^  and  \p^  as  determined  from  this  preliminary 
examination  were  then  incorporated  in  the  overall  capacity  model 

V  =  <J>,  (k,  o/ p)<|>2  (v)  +  ip  (k,  a/pH  (v)V  (2.13) 

c  i,p  z,p  i,p  2>p  b 

where  the  subscript  ’  p'  refers  to  a  given  pipe-capsule  system. 

All  available  raw  data  for  the  given  pipe  were  used  for  parameter 

estimation  and  to  evaluate  the  performance  of  the  overall  capacity  model. 

The  most  successful  relations  obtained  for  each  pipe  along  with  the 

2 

variance  of  fit  a *  and  the  95%  limits  of  the  predictions  from  these 
relations  are  given  in  Table  2.6.  The  95%  confidence  limits  (+&) 
can  be  interpreted  in  either  of  the  two  ways,  viz.,  that  the  given 
correlation  generates  95%  of  the  data  within  the  accuracy  limits  of 
+£  or  one  can  be  95%  confident  that  a  prediction  made  by  a  given 
correlation  is  correct  within  +£.  This  of  course  assumes  that  the 
prediction  is  made  within  the  region  spanned  by  the  data. 

To  simplify  the  presentation  in  Table  2.6  we  define: 

K1  =  (1-k) 


k2  =  a-k2)  e2^-) 


a 


V 

' 

■ 

' 


TABLE  2.4 


MODELING  FOR  <J>  FOR  1  1/4M  COPPER  LINE 
CARRYING  CYLINDERICAL  STEEL  CAPSULES 


Model  Form 

Variance  of  fit 

Parameters 

2 

a* 

C1  C2 

ciC)(1-k2)2 

.3022 

-21.5 

C1ln(|-)(l-k2)Z 

.03953 

-24.59 

C,  In  (2)  (1-k) 

1  P 

.129 

-13.1 

C,  ln(“)  (1-k2) 

1  P 

.149 

-  7.05 

C,ln(f) (1-k) 2 

1  P 

.043 

-81.1 

C1  In fy  (1-k4) 

1  P 

.190 

-  4.04 

1  p  (1+k2) 

.114 

-12.19 

1  p  (1+k2)2 

.0486 

-68.06 

Cii„(p) 

1  p  (1+k2) 

.0423 

-40.99 

C,  In  (2-)  (1-k)  4 

1  P 

1.228 

-20.40 

C  (2=2.)  (1-k) 

1  0 

.0832 

-20.80 

C  ln(2±2)(l+k2)2 

1  a 

1.126 

-  .386 

ci<i>(w2)2 

.757 

-  .44 

C.  (2^2)  (1-k2)2 

1  a 

.0302 

-37.604 

(Continued) 


' 
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TABLE  2.4 

(Continued) 

Model  Form 

Variance  of  Fit 

2 

a* 

Parameters 

°1  C2 

ciCP)(1-k2)2 

.302 

-21.5 

C,  (r~)  (1-k2)  2  (1+k) 

1  a 

.0285 

-20.72 

C,  (~^-)  (1-k2) 2 (1+1) 2 

1  a 

.02789 

-11.404 

c  (l-k2)  (1-k4) 

.02785 

-22.56996 

c,  c^)C2d-k2)a-kA) 

l  a 

.02689 

-21.9835  .9590 

cX)  2(l-k2)(l-k4) 

1  p 

.0495 

3.7  1.25 

C1(2^)(l-k2)(l-k4) 

.1325 

-  8.05 

1.25 
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TABLE  2.5 

MODELING  OF  ij>  FOR  1  1/4"  COPPER  LINE 
CARRYING  CYLINDERICAL  STEEL  CAPSULES 


Variance  of  Fit  Parameters 


Model  Form 

2 

a* 

di 

d2 

di(t)(p)d2(1~k) 

.00044 

.996 

.4540 

vK)d2(1'k2) 

.00047 

.996 

.2458 

1  +  d  (i  -  l)(2-)d2(1‘k  ) 

I  k  p 

.000317 

.9977 

1.199 

1  +  d1  (f  -  1)  (-)d2('1_k) 

lk  p 

.000319 

.99919 

2.1605 

1  +  d  (i  -  l)(-)d2(k  ‘  1) 

1  k  p 

1  +  d.  (l-k)(-)d2(k  "  1) 

1  P 

.000328 

1.0033 

1.7216 

.000310 

1.19945 

1.86473 

1  +  d  (l-k2)(-)d2(k  "  1) 

1  P 

.00037 

.65 

1.94 

1  +  d1(l-k)(i)d2(1'k) 

.000353 

1.196 

2.307 

1  +  (|  -  l)(j)dl(1‘k) 

.000325 

1.9937 

,£.d2(l-k) 

14-1  P 

nn4Sft 

1.098 

.694 

1  k<l-k)  + 

+  kZ] 

d  (-)d2(k  "  1) 

.612 

1  P 

0041 

1.096 

i}  ka-k)  + 1|  d-k)2  + 

rs  • 

k  J 

Table  2.5  Continued 


' 
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TABLE  2.5 

(Continued) 

1  +  d1(l-k)  (^-£-)d2(1  k) 

.00067 

1.611 

.6244 

1  +  d.  (1-k)  (^^->d2 (1-k) 

1  a 

.000955 

1.7195 

.6744 

1  +  d1(l-k2) (£zP)d2(1"k) 

.000913 

.875 

.603 

1  +  d,  (1-k)  (2n£.)d2  ^1_k  > 

1  P 

.00066 

1.612 

.3459 
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-k2'*  n-vS  ( l-a.  P 


K3  =  (l-OU-lc 


K4  =  (l-k  ) (l-k)  (-La-£-l  ■) 

a 


3/2 


K5  =  (l-k") 


K6  = 


0 


2. 3 (l-k) 


SI  =  | a-0. 8p I  + 


S2  =  |-  -  k|  +  — 
P  i  a 


|  a-0.8p 

1 


-  -k 

P  1 


where 


a-p I  =  absolute  value  of  (a-p) 


2.5  Towards  a  General  Correlation 

A  general  capacity  correlation  that  is  valid  over  a  wide  range 
of  viscosities  and  pipe  diameters  would  be  of  utmost  practical 
importance.  Although  such  a  correlation  could  not  be  obtained,  a  general 
model  form  has  been  evolved  in  the  following  few  pages. 

2.5.1  Dimensional  Analysis 


The  relation 


V  =  f(Vu,  D,  d,  a,  p,  v,  L  ,  n,  g) 
c  b  c 


may  be  approximated  by 


V  Xli  X2i  X3i  X4i  X5i  _  X6i  X7i  X8i  X9i  n  .. 
V  =  T  D  d  a  p  L  v  g  n  (2.14) 

r  h  C 

C  i=l 


where  n  is  a  suitable  integer  number  and  ^2^,  •••  are  arbitrary 

Now,  each  term  on  the  right  hand  side  of  (2.14)  must  have 


constants. 


* 


> 

/  . 


. 
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the  dimensions  of  V  .  Dropping  the  second  subscript  * i *  for  convenience, 

xl  x9  x9  x/  xc  xa  x7  Xq  xq 

[VI  =  [V.  ]  [D]  [d]  [a]  [p]  [L  ]  [v]  7  [g]  8  [n]  9  (2.15) 


where  [Z]  represents  the  dimensions  of  Z.  Thus 


[Vb]  - 

[V  ]  =  LT  1 
c 

[Lc] 

=  [d]  -  [D]  -  L 

[v]  = 

T'V 

[a] 

=  [p]  =  ML"3 

[g]  = 

-2 

LT 

M 

=  dimensionless 

The  solution  of  (2.15)  obtained  by  balancing  the  three  fundamental 
dimensions  of  M,  L  and  T  is 


x5  =  -x4 


X1  =  1  -  x7  "  2x8 


x2  =  ~x3  -  x6  -  x7  +  x8 


Therefore  (2.15)  can  be  rewritten  as 

V,  D  "X7  V3  "X8  ,  X3L  X6  X4  X9 

vc“V-T->  (Fg>  $  - 


(2.16) 


or  generally 


n  V  D  ali  W,  a2i  a3i  L  a4i  a5i  d  a6i 

\  =  I  V-f)  <rr>  $  ^  '  $ 

i=l  D  g 


(2.17) 


where  a^,  a2±,  ...  a^±  are  determined  by  x-^. 


*  *  *  x9  i  * 


,  «  I  . 
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2*5.2  Incorporating  Empirical  Observations 


Relation  (2.17)  has  been  obtained  from  purely  dimensional 
considerations.  Empirical  information  is  now  needed  to  simplify  this 
large  equation.  We  know  from  the  results  of  section  2.4  that  the 
successful  capacity  correlation  should  be  of  the  form 


vc  =  a/p,  Lc/D»  n) 

+  b2^2 (k»  Lc/D» 

+  b3^3^k,  a/p,  Lc/D,  n)Wb 
°r  vc  ■  bl  -5^5  +  b2*2Vb  +  b3*3vVb 

V 


(2.18) 


(2.19) 


Furthermore 

1)  The  second  term  may  also  be  a  function  of  viscosity.  We  do  not  have 
enough  viscosity  variations  in  the  data  to  confirm  or  reject  this 
hypothesis. 

2)  For  the  1/2"  line,  with  high  v,  the  first  two  terms  of  (2.18) 
practically  vanish. 

3)  For  other  lines  of  higher  D  and  lower  v,  the  last  term  is  not 
significant . 

A  form  of  (2.17)  that  is  consistent  with  all  the  empirical 


observations  noted  above  is 


. 

* 


' 
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where 


or 


or 


where 


VC  ■  WvVai  ♦!  +  W>*2  *2V 

D  b  V 


+  b3(-l“)1/2  *3vb 

D  g 


=  b2D-  -“1  g*  “1  +1  +  b^/2  *  V 

V 


D2  3al  e1_al  M2al-1 


<t>. 


+  b3  3/2  1/2  vVb 
D  g 


(2.20) 


>  0  and  a^  <  0.5  and 


[$±  =  4>± (k ,  a/p,  L  /D,  n)] 


Furthermore,  we  know  that 


<f>x  =  K3C1  +  C1S1k)f1(n)  (2.21) 

^  =  (K4)(l  +  C1s1k)f1(n) 

<f>2  =  1  +  C2(l-k)(^)C3(X"k)f2(n)  (2.22) 

<t>2  =  l  +  c4(i-k2)(^-)f2(n) 

<f>3  =  f3(n)  (2.23) 

1 ,  C2,  C3  and  are  positive  constants  and  f^,  f2,  are 


suitable  functional  forms. 
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2.5.3  Flow  Regimes 

What  is  the  dimensionless  nuirtber  that  determines  the  nature  of 

relationship  between  V  and  V,  ?  We  have  at  least  two  different  flow 

c  b 

regimes  corresponding  to  the  relations  (2.12a)  and  (2.12b).  One  physical 
quantity  that  apparently  does  affect  the  relationship  is  the  viscosity 
of  the  carrying  fluid. 

To  explore  the  problem  further,  the  following  empirical  analysis 
was  carried  out.  From  the  available  data  for  the  4"  steel  line,  certain 
data  points  having  the  following  characteristics  were  selected 

v  >  19  c.s. 

V,  <  1.0  ft /sec. 
b 

These  conditions  ensured  that  the  apparently  defined  Reynolds 

numbers  V^D/v  for  the  selected  data  were  well  within  the  range  of 

Reynolds  numbers  for  the  1/2"  line  data.  These  selected  points  did  not 

,  2 

fit  the  best  correlation  for  1/2"  line  very  well  (a^  =  0.0421  as 
2 

against  =  0.00133  for  the  entire  l/2"line) .  This  suggests  that  the 
Reynolds  number  does  not  determine  the  nature  of  tha  capacity  correlation. 
Furthermore,  these  selected  30  points  correlated  quite  well  with  the 
most  successful  model  for  the  4"  steel  line.  Table  2.7  compares  the 
results  of  parameter  estimation  for  the  small  subset  of  data  containing 
only  the  low  velocity  and  high  viscosity  points  against  those  obtained 
from  all  the  data. 

Considering  that  the  low  velocity  data  cover  an  extremely  small 
range,  the  agreement  of  the  parameter  estimates  in  Table  2.7  is  quite 
satisfactory.  Thus,  it  may  be  concluded  that  the  model  under  consideration 
is  a  reasonable  approximation  to  the  true  relationship.  Moreover,  it 


- 

>  ■  V\ 
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TABLE  2.7 

MODEL  PERFORMANCE  FOR  LOW  VELOCITY  DATA 
FOR  4"  STEEL  LINE 


V 

c 


K5 

0.25  2 


K5  S2 
0.25 


Low  Velocity  30  Data  Points 

All  448  Data  Points 

2 

a* 

0.004773 

0.0117 

bi 

-19.71  +  2. 749+ 

-14.21  +  0.6337+ 

b2 

-3.964  +  0.5844 

-3.208  +  0.1937 

b_ 

0.2997  +  0.0961 

0.08448  +  .00484 

3 

95%  confidence  limits. 


A3 


also  implies  that  the  nature  of  relationship  between  and  is 

independent  of  V,  .  Hence  the  dimensionless  number  that  would  determine 

b 

such  a  relationship  should  not  contain  V^.  A  dimensionless  number 

that  contains  v  and  D  but  not  V,  is 

b 

S  =  (^f)  <2-24> 

V 

* 

This  number  has  been  defined  as  the  Galileo  number  but  we  will  call  it  the 

stationary  number  sincd’  it  is  independent  of  the  fluid  velocity.  We  have 
seen  the  introduction  of  the  stationary  number  to  make  the  terms  of 
(2.20)  consistent  with  the  experimental  observations.  As  an  example, 
for  the  1/2"  line 

v  2  40  cs. 

D  =  1.3  cm 

2 

g  =  981  cm/sec 
thus  S  -  104. 

Table  2.8  gives  the  stationary  number  for  the  various  pipelines 
whose  data  have  been  used  in  this  work.  We  know  that  for  1/2  line 
with  S  -  104,  the  capacity  correlation  is  of  the  type 

Vc  *  bl  +  VVb 

and  for  other  lines  with  S  >  3x  107  the  correlations  are  of  the  type 


where  b£  is  not  a  function  of  v.  However,  there  is  too  wide  a  gap 
between  the  two  values  which  can  only  be  investigated  through  further 
experimentation. 

c.f.  Chem.  Eng.  Progress  J55,  No.  9,  p.  57  (Sept.  1959). 


4 
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TABLE  2.8 

STATIONARY  NUMBERS  FOR  VARIOUS  LINES 


Pipe  Line 


S 


1/2M  copper 
1  1/4M  copper 
4”  Steel 
4"  Plastic 


1  x  10 
5  x  10 
3  x  10 
1  x  10 


4 

8 

7 

10 

11 


10"  steel 


1  x  10 
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A  note  of  caution  seems  appropriate  at  this  stage  regarding 

the  role  of  the  stationary  number  as  a  determinant  of  capacity 

correlations.  An  equally  important  and  perhaps  even  dominant  factor 

is  the  capsule/pipe  friction  factor  n.  Again,  the  only  way  to  identify 

the  effects  of  each  factor  (stationary  number  on  the  one  hand  and 

ri  on  the  other)  separately  is  through  experimental  data. 

Furthermore,  the  stationary  number  does  not  characterize  the 

flow  mechanism  of  the  capsule  system.  All  that  is  implied  in  the 

foregoing  section  is  that  and  follow  the  same  relationship 

and  both  increase  similarly  with  the  changing  mechanism  of  flow  as 

Vfo  increases.  Thus,  S  cannot  be  used  to  determine  the  relationship 

between  PGC  and  V,  . 

b 

2.6  Pressure  Gradients 


The  same  model  building  techniques  were  used  to  try  and  obtain 
pressure  drop  correlations.  However,  no  correlation  could  be 
established  that  was  general  enough  to  be  of  any  practical  utility. 
Part  of  the  difficulty  resulted  from  the  relatively  more  complex 
nature  of  the  pressure  drop  model,  but  most  of  the  difficulty  arose 
from  the  change  in  experimental  conditions  from  run  to  run,  the 
explanation  of  which  will  be  suggested  later. 

An  examination  of  several  typical  plots  of  PGC  versus  such 
as  in  Figures  2.2  and  2.3,  suggests  very  strongly  that  some  other 
uncontrolled  parameter  has  a  much  greater  effect  on  PGC  than  . 

That  other  parameter  could  be  the  surface  conditions  of  the  capsule 
and  of  the  pipe.  This  subject  will  be  further  explored  in  the  next 


chapter . 


- 

' 

■ 
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PG  C  (psi) 
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Vb  (ft/sec) 


Figure  2.2  PGC  versus  Vb  for  run  436 

I  1/4"  copper  line  carrying  steel 
capsules  in  water.  Research 
Council  of  Alberta. 


'  >■ 


' 
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Vb  (ft /sec ) 


Figure  2.3  P6C  versus  Vb  for  run  #  450, 

I  1/4"  copper  line  carrying  steel 
capsules  in  water.  Research 


Council  of  Alberta  . 


' 
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Fating  the  data  of  a  few  experimental  runs  to  the  following 
simple  smoothing  models  provided  an  estimate  of  measurement 


variance  for  the  pressure  gradients. 


PGC  =  b ’  +  b ’ V 
1  Zb 


(2.25) 


PGC  =  b!.'  +  b"V,  +  b"V,  2 
1  zb  J  b 


(2.26) 


Table  2.9  gives  the  lower  of  the  two  variance  estimates  for  1  1/4" 
copper  line  carrying  steel  capsules.  It  can  be  seen  that  measurements 
are  becoming  more  and  more  uncertain  as  a  increases  and  k  decreases. 
Compared  to  the  small  range  of  PGC,  the  errors  are  large.  These 
errors  may  be  tolerable  for  an  individual  run  but  when  the  estimates 
of  the  parameters  from  these  runs  have  to  be  modeled  in  terms  of 
k  and  a,  one  might  end  up  with  a  totally  inadequate  or  even  misleading 
model.  For  most  of  the  pipe  lines,  even  a  qualitative  relationship 
could  not  be  obtained  between  the  shape  of  PGC  versus  and  k  and  a. 

This  can  be  seen  from  the  example  plots  in  Figures  2.2  and  2.3  for 
1  1/4"  copper  line. 

Ellis  and  Kruyer  (8)  have  attempted  to  correlate  the  pressure 
gradient  and  the  liquid  velocity  through  an  equivalent  velocity 
defined  as 


(2.27) 


V  =  V  -  A  V 
eq  an 


c 


where  A  is  an  arbitrary  constant  parameter  and  V  is  the  velocity 

cin 

in  the  annulus.  Consequently  Reynolds  number  is 


Re 


(D-d)(Van  -  A  Vc) 


(2.28) 


V 


i 

-  ' 


s 


. 


■ 
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TABLE  2.9 

PRESSURE  DROP  MEASUREMENT  VARIANCE 
FOR  1  1/4"  COPPER  LINE 
CARRYING  STEEL  CAPSULES 


Run  No. 

k 

a 

408 

.9400 

1.031 

418 

.8895 

1.125 

428 

.8390 

1.25 

436 

.8390 

1.5 

442 

.8895 

2 

448 

.9400 

3 

450 

.8895 

3 

454 

.7916 

3 

456 

.9400 

5 

458 

.8895 

5 

Range  of 

PGC 

Variance 

2 

a. 

* 

95% 

Limits 

.0  -  .03 

9.0  x  10”7 

+  .002 

-6 

.02  -  .035 

1.1  x  10 

+  .002 

.045  -  .06 

4.5  x  10"6 

+  .004 

.07  -  .085 

6.2  x  10~6 

+  .005 

.16  -  .20 

9  x  10~6 

+  .006 

-4 

.32  -  .45 

1.3  x  10 

+  .023 

.2  -  .45 

5.7  x  10"3 

±  -15 

.17  -  .30 

4.8  x  10"3 

+  -14 

.5  -  1.0 

1.9  x  10-2 

+  .28 

.5  -  1.0 

1.3  x  10~2 

+  .23 

'■ 

‘ 
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A  semi— empirical  justification  may  be  given  here  for  this 
kind  of  equivalent  velocity.  Capsule  flow  in  a  pipe  involves  two 
relative  velocities  that  contribute  to  the  viscous  drag! 

a)  annular  velocity  of  the  fluid  relative  to  the  pipe  wall: 

V 

an 

b)  annular  velocity  of  the  fluid  relative  to  capsule  wall: 

Van  **  V  SOme  kind  of  an  average  V  is  implied  here. 

The  hydraulic  radius  is  (D-d)/4.  Assuming  two  fields  are  linear 
i.e.,  the  two  velocities  can  be  linearly  superimposed  to  determine  the 
net  viscous  drag 


V 

eq 


(l-A)V  +  A(V  -V  ) 
an  an  c 


=  V  -  A  V 
an  c 


where  A  is  simply  a  weighting  factor  which  may  or  may  not  be  a  function 
of  k,  o,  D,  p ,  v. 

It  may  be  noted  that  besides  making  the  assumption  of  linear 

s uperimposition ,  Ellis  and  Kruyer  are  assuming  A  to  be  an  absolute 

constant  independent  of  k,  a,  p,  etc.  and  that  there  is  no  form  drag. 

The  last  assumption  is  not  unreasonable  in  the  light  of  our  results 

2 

of  the  preceeding  section.  The  absence  of  terms  containing  in 
the  capacity  relations  implies  no  significant  form  drag. 

However,  this  and  other  modeling  techniques  can  be  usefully 
employed  only  when  more  reliable  data  under  controlled  conditions  can 
be  generated.  In  the  meanwhile,  we  turn  our  attention  to  examining 
the  capsule  flow  mechanism  in  somewhat  greater  detail  in  the  next  chapter. 


.. 


' 


' 
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3.  A  THEORETICAL  MODEL  FOR  THE 
LAMINAR  PIPE  FLOW  OF  CAPSULES 


3.1  Scope 

The  mechanism  of  capsule  flow  in  a  pipe  is  extremely  complicated 
due  to  the  large  number  of  variables  affecting  the  flow.  What  has 
become  obvious  now  is  that  no  simple  mechanism  can  adequately  explain 
the  relationship  among  the  variables.  In  spite  of  over  a  decade  of 
effort  in  this  area,  no  satisfactory  solution  is  available  for  designing 
the  capsule  pipelines.  Limited  success,  such  as  in  Chapter  2,  has 
been  met  in  obtaining  capacity  correlations  but  the  problem  of  predicting 
pressure  drop  in  the  capsule  system  remains  largely  unsolved.  An  exam¬ 
ination  of  pressure  drops  suggests  the  presence  of  variables  that  have 
not  been  taken  into  account. 

What  is  required  is  an  examination  of  the  physics  of  the  capsule 
flow  problem.  A  significant  body  of  results  are  available  to  test  any 
theoretical  hypotheses.  For  example,  it  is  known  that  the  capsule 
normally  travels  at  a  slight  angle  to  the  pipe  axis  with  nose-up  position 
and  that  only  a  very  small  clearance  (not  usually  observable)  exists 
between  the  capsule  and  the  pipe  bottom,  even  for  capsules  with  density 
nearly  equal  to  that  of  the  liquid.  The  effect  of  capsule  density  is 
quite  pronounced  on  the  pressure  gradient  but  that  of  the  length  is 
negligible  over  the  range  of  experimentation  (length/diameter  ratio 
more  than  6).  Other  observations  are  noted  in  section  1.1. 


s 

p\  * 


' 
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3 . 2  A  Preliminary  Physical  Model 

One  can  visualize  a  capsule  being  conveyed  with  a  velocity  V 

c 

by  a  liquid  with  a  parabolic  velocity  profile.  Secondary  flows  are 
generated  from  the  requirements  of  continuity  and  are  illustrated  in 
Figure  3.1.  The  situation  in  Figure  3.1(b)  is  typical  of  conditions 
encountered  in  the  range  of  practical  velocities. 

However,  it  may  be  noted  that  the  forces  due  to  these  secondary 
flows  act  in  opposite  direction  at  the  two  ends  and  thus  the  net  force 
upward  or  downward  will  be  negligible,  although  there  may  be  an  impor¬ 
tant  contribution  to  the  torque  balance  on  the  capsule.  There  must  then 
be  other  forces  to  support  the  excess  weight  of  the  capsule.  The  analogy 
with  the  slipper  bearing  is  obvious  and  has  long  been  recognized  (5). 

In  the  slipper  bearing,  the  slipper  moving  at  a  slight  angle  to  the 
bearing  with  a  wedge  of  liquid  in  between  generates  enough  pressure  in 
the  liquid  film  to  support  its  weight.  In  the  capsule  flow,  the  capsule 
surface  could  be  visualized  as  a  series  of  slippers  of  infinitesimal 
width  moving  at  a  slight  angle  to  the  pipe  surface  (of  course  both  the 
clearance  and  the  angle  differ  with  each  infinitesimal  surface) . 

Thus  the  liquid  flow  field  could  be  considered  hopefully  to  be  the 
linear  superimposition  of  the  two  independent  flows:  (a)  Couette  flow  due 
to  the  liquid  being  carried  along  the  moving  capsule,  (b)  the  pressure 
flow  of  liquid  flowing  through  the  annulus  formed  between  the 
stationary  capsule  and  pipe  surface. 


' 


, 

' 
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SSSSSQ3SSS 

rj  ^  ^ 

X 

Xr 

Vc  * 

(a)  Vc  is  small  compared  to  liquid  bulk 
velocity  Vb. 


b 


(b)  Vc  is  of  the  order  of  Vb. 


Figure  3.1  Secondary  flows  around  the 


capsule . 
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3. 3  Couette  Flow 

3*3.1  Plane  Slipper  Bearing 

Prandtl  (25)  derived  the  motion  of  a  parallel  slipper  on  a  guide 
by  making  the  following  simplifying  assumptions  (Figure  3.2). 


(a)  p 


= 

dx 


»  the  pressure  gradient  in  the  axial  x  direction 


is  independent  of  y  on  account  of  the  small  clearance  (h) 
between  the  slipper  and  the  guide. 

(b)  Inertia  terms  are  negligible. 

(c)  Conditions  change  slowly  along  the  x-direction, 

.  92u  32u 

i.e.  -^2"  can  be  neglected  compared  to 

where  u  is  the  liquid  velocity. 

(d)  No  leakage  of  fluid  from  the  sides.  This  amounts  to 
assuming  negligible  secondary  flows  in  the  capsule  flow. 


Consider  a  stationary  slipper  of  unit  width  and  a  guide  of 
length  £  moving  in  the  opposite  direction  with  a  velocity 
(Figure  3.2).  With  the  above  assumptions,  the  equation  of  motion  for 
the  liquid  in  the  interspace  reduces  to 


92u  _  d£ 
U  3y2  "  dx 


0.1  ) 


where  p  is  the  liquid  viscosity. 

Integrating  and  evaluating  the  constants  from  the  boundary  conditions 

u  =  2_  (y2  -  ^  (hx  -  y) 

h  being  the  clearance  at  any  point  x  . 

X 


(3*2) 


' 
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x 


stationary  slipper. 
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Q  = 


fhx 

u  d  y 
o 


V  h  p'  h  3 

C  X  _  X 

2  12  vi 


This  gives 


12  y 


V 

c 

2h  2 
x 


_Q_ 

h  3 
x 


(3.3) 


(3.4) 


In  case  of  the  bearing,  the  pressure  at  x  =  0  and  x  =  £ 
must  equal  the  ambient  pressure  p^  and  thus  the  pressure  must 
rise  and  go  through  maximum  as  x  increases  from  0  to  £  . 

This  has  to  be  the  case  to  balance  the  load  carried  by  the  slipper. 
The  same  holds  true  for  the  capsule  flow  ignoring  for  the  present 
the  pressure  gradient  PGC  due  to  the  pressure  flow.  From  the 
geometry  of  Figure  3.2. 


h  - 

X 

(a- 

-x)  6 

(3 

5) 

where 

5  is 

the 

wedge 

angle. 

The  excess  pressure  at 

any  x 

is 

given 

by 

•x 

P  “  P 
x  o 

• 

p' 

o 

dx 

(3 

6) 

Using  (3.5)  in  (3.6) 


P 


x 


6px 

a62 (a-x) 


V 

c 


Q(2a-x) 

a  6 (a-x) 


(3.7) 


but  p  =  p  ,  the  ambient  pressure  at  x  =  £  ,  substituting 


' 


m 


' 
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in  (3.7)  gives 


V  a  6(a-£) 
c 


(3.8) 


(2a-£) 


and 


6PV  x(£-x) 
c 


(3.9) 


Px  Po  h  2  (2a-£) 
x 


Figure  3.3  gives  an  illustrative  plot  of  pressure  over 
the  length  of  the  guide  as  well  as  the  form  of  the  velocity 
profiles  in  the  interspace. 

Prandtl  obtained  simplified  approximate  expressions  for  the 
mean  pressure  and  friction  on  the  slipper  by  assuming  parabolic 
pressure  distribution.  We  will  proceed  further  by  using  the 
expression  in  (3.9). 

The  total  vertical  force  per  unit  width  on  the  guide  and 
the  slipper  is  given  by 


w 


(3.10) 


Substituting  for  (p  -  p  )  from  (3.9),  the  above  integral  can 

X  o 

be  reduced  to 


6yV 

c 


{ (2a-£)  £n  (-^r)  -  21} 
a-£ 


(3.11) 


If  a  >>  £,  ( %! a- £)  is  small,  then 


i 

a-£ 


(3.12) 


* 


I. 


' 
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Figure  3.3  Pressure  and  velocity  profiles  for 


slipper  bearing. 
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Substituting  in  (3.11) 


6yV  a2 

c _ 

62(2a-jO(a-*) 

3J12uVc 


where  h  =  h„ 

m  X./2 


clearance  at  the  midpoint. 


(3.13) 


Force  of  friction  on  the  upper  surface  (capsule)  will  be 
positive  and  given  by 


f3ul 


x 


Substituting  for  p’  from  (3.4) 


(3.14) 


T 

C 


A  _  JL 

2h  h  2 
x  x 


and  using  (3.8)  for  Q 


6yV  a  6  (a-£)  2yVc 
Tc  =  (2a- il)  hx2  h7 


(3.15) 


(3.16) 


The  total  drag  on  a  unit  width  of  the  guide  (capsule)  is 


f 

c 


t  dx 
c 


0 


(3.17) 


. 
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Using  the  relation  in  (3.5)  for  h  and  solving  the  above  integral 

X 


f  = 


3yV  £  y  V  (2a-£) 
— T  '  +  - r - 


£n 


m 


m 


a-£ 


(3.18) 


Again,  if  a  >>  SL  (3.18)  simplifies  to 


3yV  £  2yV  £ 

f  ~  - 9 - 9_ 

C  h  h 

m  o 


(3.19) 


Since  h^  <  hQ  ,  f^  will  always  be  positive. 

It  may  be  interesting  to  note  that  Prandtl  obtained 

yVc/hm  as  an  expression  for  the  mean  shear  stress  assuming  a 

trapezoidal  form  for  the  shear  stress  distribution.  Equation  (3.19) 

also  reduces  to  Prandtl' s  form  if  the  assumption  of  a  >>  £  is 

carried  to  its  logical  conclusion,  i.e.  h  =  h  .  In  general, 

mo 

(3.19)  predicts  somewhat  higher  drag. 

The  assumption  a  >>  £  implies  that  the  clearance  varies  only 
gradually  along  the  length  of  the  capsule  or  slipper  bearing.  This 
seems  reasonable  in  the  light  of  actual  experimental  observations  of 
the  capsule  flow  over  wide  range  of  liquid  viscosities,  capsule 
densities,  diameter  ratios  and  liquid  flow  rates.  We  shall  also  use 
this  assumption  at  a  subsequent  stage  of  this  analysis. 

From  (3..  13),  it  is  apparent  that  the  vertical  thrust  due  to  the 
pressure  distribution  in  the  clearance  is  approximately  inversely  pro¬ 
portional  to  the  square  of  the  clearance.  This  implies  that  the 
bearing  effect  in  the  capsule  mechanism  is  only  significant  over  a 
small  lower  portion.  Thus  all  the  qualitative  features  of  this  flow 
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mechanism  can  be  inferred  from  a  parallel  slipper  bearing.  Furthermore, 
since  the  flow  in  this  small  lower  section  of  the  capsule  system  will 
generally  be  laminar  regardless  of  the  nature  of  flow  in  the  upper 
layers,  these  inferences  hold  for  all  flow  regimes. 

For  example,  the  effect  of  roughness  of  both  the  capsule  and  pipe 
surfaces  on  the  pressure  drop  in  the  capsule  flow  would  be  to  cause  a 
series  of  sudden  bumping  of  the  rear  end  of  the  capsule,  i.e.  where  the 
clearance  is  least.  This  means  h  will  be  increased  at  each  bump. 

X/ 

From  relation  (3.13),  it  is  clear  that  to  support  the  same  excess  weight 

either  V  should  increase  or  the  mean  clearance  should  decrease.  The 
c 

former  would  require  large  accelerating  forces  from  the  liquid  flow 
field  to  act  on  the  capsule  and  consequently  does  not  seem  likely. 
Furthermore,  versus  plots  for  diverse  capsule  flow  conditions 

have  been  found  to  be  relatively  smooth.  Thus  the  mean  clearance 
between  the  capsule  and  the  pipe  bottom  must  decrease  resulting  in 
increasing  drag  on  the  capsule  surface  and  hence  an  increase  in  PGC 
(in  equation  (3.19)  increased  also  leads  to  similar  results). 

Theoretically,  the  two  surfaces  of  capsule  and  pipe  should  never 
come  in  contact,  for  that  would  mean  zero  clearance  which  would  generate 
very  large  pressures  that  cannot  be  balanced  by  the  excess  weight.  This 
would  be  true  if  there  were  no  secondary  flows  (side  leakages),  absolutely 
smooth  surfaces,  completely  laminar  flow  in  the  entire  annulus  and  also 
no  end  effects.  To  the  extent  the  actual  conditions  differ  from  the 
ideal,  the  motion  of  the  capsule  would  be  characterized  by  random 
deviations  from  its  most  likely  path  with  some  deviations  large  enough 
to  cause  bumping  of  the  two  surfaces.  The  secondary  flows,  as  long  as 
they  are  uniform,  cause  the  mean  clearance  to  be  consistently  less  than 
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that  predicted  by  the  no-leakage  slipper  bearing  analysis.  The  roughness 
of  the  two  surfaces,  particularly  for  low  viscosity  liquids,  aided  by 
turbulent  liquid  eddies,  may  cause  vertical  and  angular  swings  occasio¬ 
nally  large  enough  to  break  the  liquid  film  and  result  in  dragging 
motion.  In  this  sense,  both  the  clearance  and  angle  of  capsule  inclina¬ 
tion  may  be  regarded  as  random  or  more  appropriately  stochastic  variables. 

When  a  capsule  is  run  repeatedly  to  and  fro  on  a  pipe  surface 
without  rotation  (so  that  the  same  two  surfaces  face  each  other) ,  in  the 
light  of  the  above  discussion,  it  is  presumed  that  the  surfaces  change 
with  time  both  in  the  degree  and  in  the  type  of  roughness.  Large  pro¬ 
truding  humps  are  gradually  shaved  or  worn  off  and  troughs  filled  with 
metal  flow  from  adjoining  high  areas.  On  the  other  hand,  very  smooth 
surfaces  might  be  roughened  to  some  extent  by  the  nonuniform  wear  and 
hammering  effect  of  collisions.*  One  would  also  expect  the  tail  end  of 
the  capsule  to  be  worn  off  more  than  the  nose  end.  It  is  logical  to 
conclude  that  there  is  an  equilibrium  or  ultimate  roughness  which  is  a 
function  of  the  materials  of  the  two  surfaces  and  possibly  a  weak  function 
of  the  capsule  velocity.  However,  it  may  take  different  lengths  of 
running  times  to  achieve  an  ultimate  roughness,  say  only  a  few  hours  for 
a  capsule  train  with  a  fast  moving  low  viscosity  liquid  carrier  to  an 
almost  infinite  time  for  a  slow  moving  high  viscosity  liquid. 

The  effect  of  the  surface  roughness  will  be  negligible  for  viscous 
laminar  flow  due  to  the  larger  clearance,  and  the  absence  of  turbulent 
disturbances  which  could  cause  a  wave-like  motion  of  the  capsule. 


This  effect  was  actually  obseryed  at  the  Research  Council  of  Alberta 
by  J.  Kruyer. 
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For  low  viscosity  turbulent  flows,  end  effects  and  turbulent  eddies 
both  combine  to  cause  wide  fluctuations  in  the  pressure  drop.  Further¬ 
more,  since  the  mean  thickness  and  angle  of  inclination  of  the  capsule 
are  functions  of  the  surface  roughness,  the  pressure  gradient  would  also 
be  a  function  of  the  surface  roughness. 

Another  interesting  empirical  observation  is  regarding  the  effect 
of  mounting  collars  on  the  front  half  of  the  capsule.  In  the  experiments 
at  the  Research  Council  of  Alberta,  such  collars  have  been  observed  to 
drastically  reduce  the  pressure  drop  —  by  as  much  as  80%  (10) .  The  slipper 
bearing  model  readily  explains  the  observed  phenomenon.  A  collar  in  the 
front  half  of  the  capsule  would  result  in  a  greater  angle  both  due  to  the 
physical  presence  of  the  collar  and  the  increased  secondary  flows  at  the 
nose  end  because  of  the  additional  liquid  being  swept  forward  from  the 
slow  moving  liquid  layers  at  the  bottom.  The  increased  angle  of  inclina¬ 
tion  has  to  be  accompanied  by  an  increase  in  the  mean  clearance  to  support 
a  given  weight  of  the  capsule.  The  pressure  gradient,  being  approximately 
inversely  proportional  to  the  mean  clearance,  is  reduced  as  the  mean 
clearance  increases  with  the  incorporation  of  the  collar. 

Before  further  conclusions  can  be  drawn  from  this  model,  we  need 
to  extend  this  model  to  account  for  the  circular  shape  of  the  capsule 
and  also  to  include  the  pressure  flow. 

3.3.2  Capsule  Couette  Flow  Model 

Consider  now  a  capsule  being  pulled  by  a  force  F  with  a  velocity 
V  in  a  pipe  full  of  liquid  where  the  only  flow  is  due  to  the  motion  of 
the  capsule.  Any  two  facing  surfaces  on  the  capsule  and  pipe  of  infini¬ 
tesimal  width  along  the  length  may  be  considered  as  a  slipper  bearing  of 
given  "a"  and  "6"  .  Let  0  be  the  angle  that  the  infinitesimal 
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slipper  bearing  makes  with  the  vertical  axis.  The  points  A  and  B  are 
the  centres  of  pipe  and  capsule  respectively  (Figure  3.4).  For  our 
analysis  we  shall  consider  the  point  0  ,  the  midpoint  of  these  two 
centres,  as  the  approximate  centre  for  the  annulus.  If  e  is  the  eccen¬ 
tricity  of  the  two  centres,  r  and  R  the  radii  of  the  capsule  and  the 
pipe  respectively,  it  can  be  seen  that 

OC  =  r  +  y  cos  0  (3.20) 

OP  =  R  -  j  cos  0  (3.21) 

The  clearance  at  the  mid  point 

PC  =  h  (0 )  =  OP  -  OC 

=  R  -  r  -  e  cos  0  (3.22) 

Now,  the  total  drag  Fc  on  the  capsule  is  obtained  simply 
by  integrating  the  drag  forces  on  a  series  of  infinitesimal  slippers. 

[2tt 

F  =  f  (r)  d0  (3.23) 
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Using  relation  (3.19)  for  f 
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Using  relation  (3.22)  we  show  in  Appendix  A  that 
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Figure  3.4  Endview  of  capsule  and  pipe. 
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The  capsule  flow  parameters  a  and  6  (or  e  and  e  )  are 

m  o 

not  free  to  take  any  values.  They  are  constrained  by  the 
vertical  force  balance  and  torque  balance  (horizontal  force 
balance  determines  and  PGC  relationship).  The  net  vertical 

component  of  the  pressure  force  is 

r  2tt 


W  = 


w  cos  0  (r)  d0 


(3.26) 


Using  relations  (3.13)  and  (3.22)  it  is  shown  in  Appendix  B  that 
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The  pressure  distribution,  somewhat  like  that  shown  in 
Figure  3.3  ,  generates  a  torque  in  the  vertical  plane.  The  flow 
is  symmetric  in  the  cross  direction  and  only  the  force  components 
in  the  vertical  direction  need  be  considered. 

With  reference  to  Figure  3. A,  the  vertical  component  of 
force  on  an  element  (dx)  (r  d0)  at  C  is 


d2W  =  (p  ~  P  )  (dx(r  d0)}  cos  0 
x  o 


6yV  x(£-x) 
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h2  (2a- l) 


r  cos  0  dx  d0 


(3.28) 
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The  torque  about  A  (Figure  3.2)  is 
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6yVc  x2  (£  -  x) 

h2  (2a-£) 
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r  cos  6  dx  d0 


(3.29) 


Again,  using  the  relations  (3.5)  and  (3.22)  and  assuming  that 
the  clearance  is  small  compared  to  r  and  R  ,  we  show  in 
Appendix  C  that 
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Finally,  we  calculate  the  total  Couette  flow  Qc  of  the  liquid 

due  to  the  capsule  motion,  the  liquid  flow  per  unit  width  of 

slipper  bearing  being  given  by  (3.8)  •  For  an  infinitesimal  strip 

OC  +  OP 

at  P  (Figure  3.4),  the  width  for  flow  is  approximately  ( - ^ - )d0  • 


V  a  6  (a-£) 
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(3.31 ) 


Using  relations 

dQ  = 
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(3.5),  (3.20),  (3.21),  and  (3.22) 
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Qc  = 
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It  may  be  noted  that  (3.33)  also  gives  the  flow  for  concentric 
capsule  flow  with  a  perfectly  linear  velocity  profile  in  the  annulus. 

3.4  Pressure  Flow 

Now  consider  a  stationary  capsule  with  liquid  flowing  through 
the  annulus  with  an  imposed  pressure  gradient  PGC  .  We  shall  confine 
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our  attention  to  the  case  where  the  flow  in  the  annulus  is  laminar. 
Neglecting  end  effects  and  secondary  flows  and  for  the  steady  state 
conditions , 


<PGC>  * 


(3.34) 


With  reference  to  Figure  3.4,  u  is  the  velocity  at  any  angle  0 
and  distance  ^  from  0. 

Integrating  twice 
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(3.35a) 
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(3.35b) 


Using  the  end  conditions  u=0  at  \p  =  r  +  —  cos  0 

£ 

and  also  u=0  at  ip  =  R  -  —  cos  0  to  evaluate  and  C ^  » 

equation  (3.35b)  reduces  to 
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where  r. 


R  -  ~~  cos  0 


(3.37a) 


and 
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r  +  y—  cos  0 


(3.37b) 


Again  we  have  used  e  for  e  and  consider  it  independent  of  x  . 
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The  total  pressure  flow  in  the  annulus  is  given  by 


f2TT 

0 


u  (\p  d  0)  dip 


(3.38) 


We  show  in  Appendix  D  that  approximately 


5tt  R4 
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(3.39) 


where  k  =  r/R  . 

The  results  in  (3.39)  may  be  compared  with  those  of  Tiedt 


(26). 


f  =  "  (PGC) 

P  u 


(3.40a) 


where  p  is  the  liquid  density  and  u  is  the  average  annulus 
velocity  given  by 


u  =  Q  /  ttR2  (1-k2 ) 
P 


(3.40b) 


The  Reynolds  Number  is  defined  as 

R  =  (D-d)  u 


(3 . 4  0c  ) 


This  gives 

f  Re  =  27TR-  q-^  (1-k2)  (1-k)2 


(3 . 4  Od  ) 
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Substituting  for 


(3.3  9)  gives 


f  R  = 
e 


64 

5 


=  12.8 


(3.40e) 


The  agreement  of  (3.40e)  with  the  results  obtained  by  Tiedt 
is  not  as  good  as  might  have  been  hoped  for.  Nevertheless,  it  is 
good  enough  for  our  purpose  considering  the  approximations  involved. 
According  to  Tiedt,  for  zero  bottom  clearance  and  k  =  1/2,  the  value 
for  f  Re  is  10.3  and  it  approaches  9.6  as  k  approaches  unity. 
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The  shear  force  on  the  capsule  surface  due  to  pressure  flow 
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From  relation  (D— 5)  Appendix  D 
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The  total  force  on  the  surface  is 
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The  moment  in,  the  vertical  plane  about  D  (Figure  3.5)  caused  by 
the  horizontal  shear  forces  on  the  surface  due  to  both  the  Couette  flow 
and  the  pressure  flow  is 


-  (f£  r  d0  +  £  r  d0)  r(l  -  cosQ) 


(3.46) 


Using  relations  (3-19)  and  (3.44)  ,  it  is  shown  in  Appendix  E 
that 
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3.5  Total  Flow 

We  now  superimpose  the  Couette  and  the  pressure  flows.  We 
assume  the  two  flow  fields  are  linear,  i.e.  their  effects  are  additive. 
From  the  requirements  of  continuity 


Total  flow  *  Couette  flow  +  Pressure  flow  +  Capsule  flow 
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Figure  3.5  An  infinitesimal  surface  on  the  capsule. 
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or 


iTR2Vb  =  Qc  +  Qp  +  irr2Vc 


(3.40) 


where  is  the  bulk  or  the  average  velocity  in  the  capsule 

flow  system.  Substituting  for  Q  and  Q  from  (3.33)  and  (3.39) 

c  p 


Vb  =  32^7  (PGC)  CL“k2>(l-k>2+  \  <1+k2>  vc  (3*49> 


The  forces  acting  on  the  capsule  must  balance.  For  example, 
if  secondary  flows  are  neglected  the  excess  weight  of  the  capsule 

=  Vertical  component  of 
the  excess  pressure  due 
to  Couette  flow. 
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When  the  clearance  is  small,  so  that  -  (R-r)  , 

the  second  term  in  the  parentheses  Ls  much  smaller  than 
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Similarly,  the  forces  in  x-direction  must  balance,  i.e. 

Pressure  force  +  viscous  drag  due  to  Couette  flow 
+  viscous  drag  due  to  pressure  flow  =  0 
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Moment  about  A  due  to  Couette  pressure  forces 
+  Moment  due  to  excess  weight  +  Moment  due  to  all 
viscous  forces  =  0 
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i.e. 
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We  now  have  five  equations  (3.49) , (3.52) ,  (3.53),  (3.54),  (3.55) 

and  five  unknowns  viz.  V  ,  PGC  ,  e  ,  e  ,  and  e  . 

c  m  3c  o 

When  the  clearance  and  the  angle  5  are  small,  it  can  be 


shown  (Appendix  F)  that 
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and  for  a  high  viscosity  liquid 
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It  may  be  noted  that  torque  balance  was  not  used  to 
arrive  at  the  above  results;  thus  the  significant  effect  of 
secondary  flows  on  the  torque  did  not  enter  our  results. 


3.6  Physical  Contact  of  Surfaces 

Frequent  bumping  of  the  capsule  into  the  pipe  surface  and  actual 

dragging  at  one  end  would  cause  an  additional  drag  force  and  an 

upward  reaction  force  on  the  capsule.  The  upward  reaction  force  may 

be  considered  as 
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and  the  horizontal  drag 


positive  in  the  x-direction,  is 
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where  n  is  the  coefficient  of  wet  friction  (i.e.  two  surfaces 
are  wetted  by  the  carrier  liquid). 

Thus  (3.50)  becomes 
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Proceeding  as  before , (3 . 61)  reduces  to 


2  ' 
e 

6  u  £  V 

c 

\ 

(R-r)2 

r (R-r)2  (a-£)  g 

• 

r 

l+4> 

f  £  \ 

m 

»  * 

[(R-r)J 

> 

2/3 


(3.62) 


Similarly  (3.51)  results  in 
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Substituting  from  (3.62) 
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Substituting  PGC  from  (3.64)  in  (3.49) 
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From  (3.59)  , 
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may  be  interpreted  as  the  fraction  of  the 


((R-r)J 

excess  weight  supported  by  the  surface  reaction  or  that  supported  by 
secondary  flows.  Thus  <J>  may  be  an  empirical  term  that  incorporates 
the  effect  of  all  simplifying  assumptions  in  our  model. 

When  <p  is  small  compared  to  unity,  i.e. 
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(3.64)  reduces  to 
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and  (3.65)  now  gives 
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3 . 7  Threshold  Velocity 

Threshold  velocity  is  the  liquid  bulk  velocity  which  just  sets 
the  capsule  moving.  Since  the  mechanism  of  flow  is  the  same,  equation 
(3.65)  should  yield  the  threshold  velocity  V,  when  V  =0  except 

D  j  U  C 

that  now  r|  would  have  to  be  replaced  by  n  ,  the  coefficient  of  static 
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wet  friction.  Thus 


b,t 


5 R2  (1-k2)  (1-k) ^  k2  g  ns 

32  y 


$ 


m 


(R-r) 


(3.69) 


3. 8  Model  Validation 

In  view  of  the  assumptions  made  in  the  derivation,  we  would  not 
expect  the  model  to  be  able  to  predict  quantitatively  the  true  relation¬ 
ship  between  the  different  physical  variables.  Hopefully,  the  model 
would  predict  the  right  trends  and  suggest  the  form  of  some  of  the 
relationships . 

We  have  selected  a  typical  run,  No.  64,  from  the  Research  Council 
of  Alberta  data,  which  involved  the  transport  of  steel  capsules  in  oil  in 
a  1/2"  diameter  copper  line.  Figures  3.6  and  3.7  give  comparisons  of  the 
observed  and  predicted  relationships  between  and  PGC  (equations  (3-56) 
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and  (3.57  ))  °n  the  one  hand  and  V,  on  the  other. 

b 

As  expected,  both  the  predicted  VL  and  PGC  for  a  given  V  are 

b  c 

lower  than  the  corresponding  observed  values.  On  examining  the  key 
assumptions  used  in  deriving  equations  (3.56)  and  (3.57)  namely,  laminar 
flow,  sufficiently  small  bottom  clearance  compared  to  R-r  ,  and  the 
absence  of  any  secondary  flows,  we  feel  that  the  least  warranted  assump¬ 
tion  is  the  last  one.  The  bottom  clearance  would  be  significantly 
decreased  by  the  presence  of  secondary  flows  (analogous  to  the  side 
leakage  in  the  slipper  bearing,  reference  25)  causing  greater  drag  on  the 
capsule.  Consequently  in  actual  pipeline  capsule  flow  with  secondary 
flows,  a  liquid  bulk  velocity  higher  than  that  predicted  by  the  model 
would  be  required  to  achieve  a  given  capsule  velocity.  Also  the  pressure 
gradient  observed  should  be  higher  than  that  predicted. 

In  addition,  surface  roughness  and  capsule  fluctuations  would 
cause  the  capsule  to  be  in  physical  contact  with  the  pipe  bottom  for  a 
certain  fraction  of  the  time.  This  would  result  in  additional  drag  on 
the  capsule  and  thus  compound  the  effect  of  secondary  flows. 

Comparing  the  results  of  this  model  (V^  versus  Vc)  with  the 
empirical  relations  obtained  in  Chapter  2  from  the  data  of  the 
Research  Council  of  Alberta,  the  two  model  forms  look  somewhat  different. 

The  empirical  results  for  1/2"  copper  line  in  laminar  flow  indicate 

is  proportional  to  whereas  equation  (3.57)  includes  two  terms,  one 

2/3 

proportional  to  V  and  the  other  proportional  to  V  ,  and  the  latter 

c  c 

term  although  smaller  than  the  first  is  significant.  However,  it  was 
found  that  the  prediction  of  equation  (3.57)  was  quite  close  to  the 
straight  line  relationship  (see  Figure  3.6  as  an  example)  found  empirically. 
Pressure  drop  data  could  not  be  relied  upon  for  a  meaningful  comparison. 
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lhe  utility  of  this  model  can  be  further  enhanced  by  an  empirical 
modification.  Incorporation  of  one  free  parameter  each  in  relations 
(3.56)  and  (3.57)  results  in  much  better  agreement  with  the  experi¬ 
mental  results.  The  predictions  for  the  modified  relations 
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are  shown  in  Figures  (3.8)  and  (3.9).  The  parameters  b^  and  b^ 
are  estimated  from  the  experimental  data  themselves  giving  the 
following  least  square  estimates 

b  =  1.714  +  0.0443  (95%  confidence  limits) 

Standard  deviation  of  fit  for  model  (3.70) 

=  0.0983  or  approximately  5%  of  the  average  prediction. 

b?  «  187.1  J- 14.08 

Standard  deviation  of  fit  for  model  (3.71) 

=  0.087  or  approximately  3%  of  the  average  prediction. 
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These  conclusions  are  at  best  preliminary  and  need  further 
work.  Some  independent  method  of  either  directly  measuring  the  surface 
roughness  or  of  inferring  these  properties  would  be  particularly 
helpful.  Measurements  of  the  actual  bottom  clearance  and  angle  can 
be  very  useful  in  establishing  the  capsule  flow  mechanism  and  the 
extent  of  the  secondary  flows. 

In  conclusion  then,  the  slipper-bearing  model  holds  promise  as 

a  framework  for  further  analysis.  Qualitatively,  it  explains  empirical 

observations  of  small  bottom  clearance,  the  slight  angle  with  nose-up 

position  for  the  moving  capsule,  the  effect  of  surface  roughness, 

the  effect  of  collars,  and  relationship  among  the  key  variables  V^, 

V  ,  and  PGC.  Further  modifications  and  development  may  also  provide 
c 

reasonable  quantitative  agreement  with  the  observed  data. 
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4.  CONCLUSIONS 


4 . 1  Summary 

The  concept  of  transporting  solids  in  the  form  of  capsules  in 
pipelines  evolved  at  the  Research  Council  of  Alberta  during  the  1950's. 
Since  then,  substantial  amount  of  research  effort  has  been  devoted  in 
this  area  with  a  view  to  gaining  sufficient  understanding  of  the 
system  for  possible  commercial  exploitation.  While  extensive  experi¬ 
mental  data  have  been  generated  at  the  Research  Council  of  Alberta  and 
elsewhere,  no  comprehensive  design  solution  is  yet  available. 

An  attempt  has  been  made  in  Chapter  2  to  statistically  analyse 
most  of  the  data  available  from  the  Research  Council  of  Alberta  to 
obtain  general  correlations  between  the  capsule  velocity  and  pressure 
gradient  on  the  one  hand  and  the  independent  variables  -  bulk  velocity, 
liquid  viscosity  and  density,  capsule  diameter,  length  and  density,  and 
the  pipe  diameter  on  the  other  hand. 

While  a  completely  general  capacity  correlation  requires  further 
experimental  and  modeling  work,  useful  correlations  have  been  obtained 
for  individual  pipelines  ranging  in  size  from  one  half  of  an  inch  to 
ten  inches.  The  effect  of  capsule/pipe  diameter  ratio,  liquid 
viscosity,  capsule /liquid  density  ratio  and  the  liquid  bulk  velocity 
could  be  modelled.  However,  the  effect  of  the  surfaces'  roughness 

and  pipe  diameter  could  not  be  separately  identified  due  to  lack  of 

( 

relevant  data. 

Little  success  has  been  obtained  in  predicting  the  pressure 
drop.  The  pressure  drop  data  not  only  had  large  variance  but  worse, 
they  were  found  to  be  generally  not  reproducible.  This  indicates  the 
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presence  of  variables  that  have  not  been  taken  into  account  in  planning 
the  experiments.  It  is  concluded  that  the  capsule/pipe  surface 
conditions  are  among  the  important  variables  that  were  not  measured. 

The  second  approach  to  the  problem  was  mechanistic  modeling  of 
the  capsule  flow.  A  preliminary  physical  model  envisaging  the  capsule 
flow  as  consisting  of  a  series  of  infinitesimal  parallel  slipper 
bearing  is  presented  in  Chapter  3.  Using  a  set  of  simplifying  assumptions 
including  one  which  neglects  secondary  flows,  capacity  and  pressure 
drop  correlations  have  been  developed  for  laminar  flow.  As  expected 
from  the  assumptions  made,  the  model  predictions  for  both  the  bulk 
velocity  and  the  pressure  gradient  for  a  given  capsule  velocity  are 
lower  than  the  observed  values.  However,  addition  of  only  one  empirical 
parameter  in  the  correlations  gives  reasonable  agreement  with  the 
observed  values.  The  approach  appears  to  be  quite  promising. 

4 . 2  Further  Work 

Considerable  amount  of  further  research  is  required  both  on 
experimental  and  modeling  fronts  before  a  complete  understanding  of 
the  basic  fluid  mechanics  of  capsule  pipeline  flow  is  possible.  Addi¬ 
tional  controlled  data,  particularly  for  pressure  drop  and  empirical 
modeling  of  data  on  the  lines  described  in  Chapter  2  are  needed.  It 
will  also  be  useful  to  incorporate  secondary  flows  in  the  theoretical 
slipper  bearing  model  developed  in  Chapter  3. 

To  obtain  a  general  capacity  correlation  applicable  to  a  range 
of  pipeline  sizes  and  capsule/pipe  surfaces,  data  should  include 
measurements  of  coefficient  of  friction.  In  addition,  a  study  of 
changes  in  capsule/pipe  surface  roughness  over  time  would  be  useful 
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for  commercial  start  up.  In  experimental  runs,  the  capsule  may  be  run 
to  and  fro  for  a  sufficiently  long  time  to  achieve  ultimate  roughness 
(as  defined  in  section  3.3.1)  before  pressure  drop,  velocity  and 
friction  data  are  recorded.  To  determine  the  viscosity  dependence  in 
the  capacity  correlation  more  precisely,  data  in  the  viscosity  range 
of  1  c.s.  to  20  c.s.  for  a  1/2"  line  or  1  c.s.  to  30  c.s.  in  a  1  .1/4" 
line  and  viscosity  varying  in  a  single  run  from  1  c.s.  to  10  c.s.  in  a 
1  1/4"  line  or  from  2  c.s.  to  30  c.s.  in  a  4"  line  would  be  helpful. 

The  model  developed  in  Chapter  3  seems  to  be  a  step  in  the  right 
direction.  It  is  expected  that  substantial  improvements  in  predictions 
would  result  if  secondary  flows  can  be  incorporated  in  the  model.  The 
secondary  flows  have  the  effect  of  reducing  the  bottom  clearance  in  the 
parallel  slipper  bearing  and  thus  give  a  higher  actual  pressure  drop  than 
that  predicted  by  the  no-leakage  bearing  model.  The  treatment  for 
incorporating  side  leakages  (analogous  to  the  secondary  flows  in  the 
capsule  system)  is  given  by  Michell  (27)  and  could  be  used  as  a 
starting  point. 

A  more  ambitious  development  of  the  slipper  bearing  model  would 
be  to  divide  the  annulus  into  two  zones  with  the  liquid  in  the  lower 
section  (of  smaller  clearance)  in  laminar  flow  and  in  turbulent  flow 
in  the  upper  section.  A  suitable  velocity  profile  may  be  assumed  for 
the  turbulent  flow.  The  angle  including  the  lower  laminar  section 
could  be  arbitrarily  assumed  to  be  a  function  of  Reynolds  number. 

While  we  are  still  a  long  way  from  the  level  of  sophistication 
achieved  in  the  fluid  flow,  substantial  progress  has  been  made  qs  a 
result  of  the  work  done  at  the  Research  Council  of  Alberta.  It  is  felt 
that  further  developments  in  the  above  direction  would  significantly 
improve  our  understanding  of  the  capsule  pipeline  flow. 
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NOTATION 

a  distance  between  the  capsule  nose  and  where  the  line  of 

capsule  surface  intersects  the  pipe  surface,  cm 
a^  dimensionless  constants  (i  =  1,  2,  3,  4,  5) 
empirical  parameters  (i  =  1,  2,  3) 
d  capsule  diameter,  cm 

f  drag  per  unit  width  of  the  guide,  dyne/cm 

2 

g  gravity,  cm/sec 

h  liquid  film  thickness,  cm 

k  diameter  ratio  d/D,  dimensionless 

2 

p  liquid  pressure,  dyne/cm 

2 

p^  ambient  pressure  dyne/cm 

2 

pT  pressure  gradient  for  the  capsule  system,  dyne/cm 

r  capsule  radius,  cm 

r^  distance  of  capsule  surface  from  the  center,  cm 

r 2  distance  of  the  pipe  inside  surface  from  the  center,  cm 

u  local  liquid  velocity,  cm/sec 

w  vertical  force  per  unit  width  of  the  slipper,  dyne/cm 

x  distance  in  axial  direction,  cm 

x.  dimensionless  constants  (i  =  1,  ...  9) 

l 

y  distance  in  vertical  direction,  cm 

C.  integration  constants  or  empirical  parameters 

(i  =  1,  2,  3,  4) 

D  pipe  diameter,  cm 

E  roughness  factor,  dimensionless 
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F 


K. 

x 

PGC 

PGF 


PR 

Q 

R 

S. 

l 

T 

V 

VR 

W 


total  drag  on  the  capsule,  dynes 

functions  of  k,  o  and  p  (1=1,  ...  6)  defined  in  section  2.4 

pressure  gradient  for  the  capsule  system,  dyne/cm^  or  psi 

pressure  gradient  when  liquid  alone  flows  at  the  rate  of 
3 

V^»  dyne/cm  or  psi 

pressure  ratio  PGC/PGF,  dimensionless 

2 

liquid  flow  rate  per  unit  of  width,  cm  /sec 
pipe  internal  radius,  cm 

function  of  a  and  p  (1*1,  2),  defined  in  section  2.4 

torque  on  the  capsule,  dyne-cm 

velocity,  cm/sec  or  ft/sec 

velocity  ratio,  dimensionless 

total  vertical  force,  dyne 


6  angle  between  the  infinitesimal  capsule  surface  and  the 

pipe  surface,  radian 

e  eccentricity  of  the  capsule  and  pipe  centres,  cm 

H  capsule  length,  cm 

3 

p  liquid  density,  g/cm 

3 

o  capsule  density,  g/cm 

p  liquid  viscosity,  poise 

v  liquid  kinematic  viscosity,  stoke  or  centistoke 

0  angle  that  the  infinitesimal  surface  makes  with  the  y-direction, 

radian 

2 

t  shear  stress  or  drag  force  per  unit  area,  dyne/cm 

ip  distance  of  a  liquid  shell  from  centre  point,  a  general 


functional  form 


. 

.  IT 

> 


> 
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<f>  a  general  functional  form 

n  coefficient  of  friction  between  capsule  and  pipe  surface, 

dimensionless 

Subscripts 

an  annulus 

b  bulk 

c  capsule 

l  the  tail  end 

m  the  mid  point 

o  the  nose  end 

P  pipe 

s  stationary 

t  threshold 

x  at  a  distance  x  from  the  nose  end 

C  capsule,  Couette  flow 

H  horizontal,  in  the  axial  direction 

P  pressure  flow 


V 


vertical 


> 


' 


' 
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APPENDIX  A 


F  =  2yV  £ 
c 


o 


LR-r-e  cos0  (R-r-e  cos0) 
ra  o 


]  r  d0 


2yV  £r 
c 

(R-r) 


-]  d0 


o  [1  - 


m 


(R-r) 


cosO]  [1  - 


(R-r) 


cosO] 


(A-l) 


But 


TT 


d0 


(1  - 


(R-r) 


cosO ) 


arctan 


(1  +  ^)tan  - 


[1  - 


(R-r)2 


] 


[1  (R-r) 2  ^ 


1/2 


o 


IT 


[1  - 


(R-r)2 


1/2 


(A-2) 


Substituting  (A-2)  in  (A-l) 


F 


2TryV  £r 

_ r 

(R-r)  1 


3 


[1 


£m  -,1/2 


[l- 


(R-r) 


(A-3) 
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APPENDIX  B 


W  =  2  J  w  cosO  (r)  dO 
o 


6 £  mV  r 
c 


cosO  d0 


o 


(R-r-e  )  cosO) (R-r-e .  cosO) 
m  £ 


6£  mV  r 
c 


[- 


dO 


“I" 


(R-r)dO 


-]  (B-l) 


o 


(R-r-e  cosO)  e. (R-r-e  cosO)(R-r-e„  cosO) 
m  £  m  £ 


Using  the  results  in  (A-2) 


W  = 


6£  mV  r 


TT 


[  " 


7T 


[1  - 


m  .jl/2 


dO 


(R-r) 2 


o  [1  - 


m 


(R-r) 


cosO] [1  - 


'£ 


(R-r) 


cosO  ] 


(B-2) 


The  expression  in  (B-2)  may  further  be  simplified  when  the  angle  6 

is  small  so  that  e  and  e.  are  nearly  equal,  i.e. 

m  £ 


[!  -  cose]  [l  - 


(R-r) 


£  Em+e£  2 

COS0]  -  [1  -  - vCOS0] 


2 (R-r) 


(B-3) 


But 


[1  - 


(t  +E.) 

m  £ 

4 (R-r) 2 


dO 


o  [1  - 


e  +e 

m  £ 

2 (R-r) 


dO 


cosO]  o  [1  - 


e  +e 

m  £ 

2 (R-r) 


cosO  ] 


7T 


(e  +  e  )  sinO 

m  x. 


e  +e 

2(R-r)  [1  -  ™  j-  cosO] 


R-r 


(B-4) 


101 


or 


(f-  +c  ) 

[1 - -  ■  -9~  ] 

4 (R-r) 


r  tt 


[1  - 


dO 


e  +en 
m  ft 

2 (R-r) 


COS0  ] 


(  e+e  n )  2  I/O 

f  i  ^  i  1/  ^ 

[1  '  4 (R-r) ^ 


(B-5) 


Hence 


/■it 


(1 - — 

y  R-r 


dO 


IT 


m  cos0)(l  -  — —  cosO) 
R-r 


d0 


o  [1  - 


e  +e 
m  ft 

2 (R-r) 


cosO] 


(e  +e  )2  3/2 

m  x,  , 

2  ■* 

4 (R-r) 


(B-6) 


Substituting  in  (B-2) 


W  = 


67Ty£  V  r 
Eft  <*■*>" 


(e  +E.)2  3/2 

n  _  m  £ _ t 

11  4 (R-r) 2  J 


e  1/2 
[1  -  —m—  0-] 


(R-r) 


(B-7) 
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APPENDIX  C 


T  = 


2tt  6 yV 


cos0  [ 


£  2, 

x  (£-x)rdx 


-]  d0 


o  x 


Substituting  for  h  from  (3.5) 

X 


(c-i) 


2  it 


o 


x^(£-x)  r 

2  2 

(a-x)  6 


dx 


~~2  [_x 

6 


-  (2a-£)  + 


a(3a-2£) 

(a-x) 


(a-£) a' 
(a-x)2 


-]  dx 


2  2 

[-  ^ - (2a-£)x  -  a(3a-2£)  £n  (a-x)  — a  ] 

2  (a-x)  o 


=  ^2  [-  T~  -  (2a"^+  3 (3a-2£)  £n  (^)  -  a£] 

6 


(C-2) 


When  the  angle  6  is  small,  i.e.  a  >>  £  ,£n( — -)  may  be  approximated 

3. 

b'y  the  first  two  terns  of  its  Taylor  series. 

2 


1  £ 


£n  (a-£)  £n(1  a-£  )  ~  (a-£)  2  (a-£)z 


(C-3) 


Substituting  in  (C-2)  and  simplifying 

x2(£-x)  r  j  r  r  „  „  .  £2  ,  a£(3a-2£)  1  a£2(3a-2£)  t 

- 2 — 2~  dx  -  ~2  [_3a£  +  j-  +  — 7--p: - 7 - ^  ] 

(a-x)  6  6  2  Ca  £)  2  (a-£) 


o 


. 


. 


103 


2  2 

_  r  r£  (3a-£)  a£  (3a-2£)n 
52  L  2 (a-£)  2(a-£)2  J 


_r_  £3 (2a-£) 

2  2 
6  2(a-£) 


(C-4) 


Substituting  in  (C-l) 


T  =  - 


2tt  ,  „  3 

6yV  r  £  cos0 
c 


o 


2  2 

2(a-£)  6 


7T 


6yV  £' 
c 


cos0  d0 


o  h£ 


,Ti 


6yV  £' 
c 


r  cos0 


o 


(R-r-e  cos0)2 

X/ 


d0 


6yV  £  r 
c 


(R-r) 


Q  ej,(R-r-e£  cos0)  G£^R-r“££  cos0) 


-]  d0 


(C— 5) 


Utilizing  the  results  in  (A-2)  and  (B-5) 


6-rryV  £  r 

T  =  - - -  [ 

e^(R-r) 


r,  1 3/ 2  r,  e£  ,1/2 

[1  “  (R-r)"2  [1  "  7wa] 


] 


67ryVc£  r 


(R-r)3[l  -  Td“ZTz]  3/3 


(R-r) 


(C— 6) 


,  '  .  • . 
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0 

P 


dO 

o 


u  ip  d  ip 


=  2 


q  d0 
P 


o 


where 


r  r. 


qP  = 


u  tp  d  ip 


(D-l) 


(D-2) 


may  be  considered  as  flow  per  unit  of  angle  around  0  (Figure  3.4) 
Substituting  for  u  from  (3.36) 


PGC 
qp  4y 


,  2  2 

rr  (r  —  r  ) 

\  ^  '“I  2 

[-  \p  H - ip  Zn  ip 


Zn 


2  2 
(r„  £n  r  -r  £n  r„) 

+  — - - 1 — - —  i|;]  d^ 


£n 


PGC  [*4  +  ^  -r-2-  (2*2  £n  *-*2) 


16y 


£n 


(r22  tn  r  -  r/  In 
+  - - -  (2 ip  ) 

Zn  — 


(D-3) 
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Taking  the  limits  and  simplifying 


2  2 
(v  —  r  ) 

PGC  ,2  2,  r  2  .  2  k  1  2  , 

qp  “  1617  (rl  -  r2)  [rl  +  r2 - r - 1 

Jin  — 
r2 


(D-4) 


The  natural  log  may  be  approximated  by  the  first  two  terms  of  its 
series  expansion  when  (r^-r^)/r^  is  small,  i.e. 


£n 


r  — r 
1  2 

!ln  [1  +  — — -] 


Vr2  1 


r2  2  r  2 
2 


rl  r2  1  rl  r2 

-i— ^  [1-4  ~J—^] 

r.  2  Tn 


(D-5) 


or 


r2 


rl  r2 


1  rrr2 -1 

[i  - 

2  r  „ 


rl  r2 


2  1  rl  r2 

[i+4  4—4 


2  r, 


- - —  +  j 

rrr2  2 


(D-6) 


Substituting  (D-6)  in  (D-4)  yields 


q  =  PGC  f  .  \  f  N  3 

P  32y  r!+r  2  ^rl  r 2 


(D-7) 


. 
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Now 


=  R  -  —  cos  0 


r2  =  r  +  -  cos  0 


R  -  m  . 
-y—  cos  0 


,  m  r. 

r  +  —  cos  0 


and 


(R-r) 


m 


Thus 


q„  = 


P  32y 


PGC  (R+r) (R-r) 3  [1 


2  3 

3cos0  +  3cos  0  -  cos  0] 


(D-8) 


0 


2 (PGC) (R2-r2) (R-r)2 

32y 


IT 


[1  -  3cos0  +  3cos  0 


cos  0]  d6 


(D-9) 


Noting  that 


cos0  d0  =  sin0 


2n  ,n  sin0  co s6  ,  1 
cos  0  d0  =  - 2 - 1  0 


(D-10) 


3„  sin0  cos  0,2  A 

cos  0  d0  =  - - - h  —  sin0 


and  that 


k-t 


(D-9)  simplifies  to 


O  =5,R  (PGC)  (1-k2  (1.k): 
‘P  32y 


(D-ll) 


\ 
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APPENDIX  E 


Substituting  f  from  (3.19)  and  t  from  (3.44)  in  (3.46)  and  using 
c  p 

the  mean  clearance  e  for  e. 

m 


2tt 


th  - 


3pV  l  2pV  i 


o 


h  h 

m  o 


-]  r^(l  -  cosO)  d0 


2tt 


(PGC)  Hr  rl 

4r„ 


[-j  (r^  +  r2  ^  -  2r22]  (l-cos0)  dO 


o 


=  2pV  Sir 
c 


■-3(l-cos0)  d0  _  2(l-cos0)  d0  , 

(R-r-e  cosO)  (R-r-e  cos0) 
m  o 


(PGC)  £  r 
8 


IT 


(R+r) (R-r-e  cos0) fl-cos0) 
m 


m 


(r  +  y~  c°s0) 


4-  2(R-r-e  cos  0) (l-cos0)  d0] 
m 


(E-l) 


Using  the  result  in  (A-l) 


TT 


7T 


(1-CQS0)  d0 

(R-r-ecos0) 


£ 


R-r-e 


(R-r-ecos0) 


■]  d0 


o 


(R-r-e)  tt 

2 

0  1/2 


(E-2) 


' 
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and 


TT 


O 


(R-r-e  cosO) (1-cosO) 


(r  +  y-  cosO) 


2(e  +R+r)  (R+r)(e  +2r) 

[2cos0 - - -  +  - - - ]  dO 


o 


m 


e  (r+~  cosO) 
m  2 


2tt(e  4-R+r)  (R+r)(e  +2r)ir 

m  m 


m 


e  1/2 

er[l+ 
m  ,  2 
4r 


(E-3) 


Also 


7T 


(R-r-e  cos0) (1-cosO)  d0 


IT 


[(R-r)  -  (R-r+e  )  cos0  +  e  cos  0 ] d 0 
m  m 


tt  (R-r  +  -f) 


(E-4) 


Substituting  (E-2) ,  (E-3)  and  (E~4)  in  (E-l)  gives 


T  =  2tt)jV  Hr' 
R  c 


3  2 


3 (R-r-e  ) 
m 


1/2 


2 (R-r-e  ) 

+ - 2 — z 


C°  ?R=FF>  '  ^o(R-r)(1-  75I7TT) 


£  “  1/ 2j 

(R-r)-' 


+ 


tt(PGC)  l  x‘ 
8 


2 (e  +R+r) (R+r) 


m 


m 


(R+r)  (e  +2r) 
m 

xl/2 

e  r(l+  -~) 
m  ,  2 
4r 


+  2  (R-r  -  y) 


(E-5) 


■ 
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APPENDIX  F 


m 


When  the  clearance  and  the  angle  6  are  small  so  that 
~  e  -  e  =  (R-r) ,  we  may  replace  e  and  e  by  e  and  e  by 


(R-r)  where  the  expression  is  not  very  sensitive  to  c  . 

m 

Thus  (3.52)  reduces  to 


[1  ~ 


m 


(R-r) 


rl  =  [ 


6y  £V 


2/3 


emr(o-p)  g(R-r) 


(F-l) 


For  example,  the  R.H.S.  is  not  sensitive  to  the  value  of 

and  may  be  replaced  by  (R-r)  as  an  approximation  but  the  L.H.S. 

is  very  sensitive  to  the  value  of  e  (e  =  R-r  would  reduce  it  to 

mm 

zero) .  The  effect  of  such  approximation  depends  upon  the  value  of 

£  ,  £  and  £„.  As  an  example,  a  bottom  mean  clearance  of 
m  o  i  ’ 

h 

[h  /(R-r)]  =  0.0001  and  h.  =  would  introduce  an  error  of  around 
m  SC  Z 

one — hundredth  of  one  percent  in  approximating  relation  (3.52)  to 
(F-l)  . 


Thus 


1 


m 

(R-r) 


7 


6m  IV 

^ (R-r) zr (a-p) 


(F-2) 


or 


e  =  (R-r) 
m 


6  m£V 


1  *-(R-r)zr(a-p)  g^ 


2/3 


(F-3) 


Similarly  (3.53)  reduces  to 


16pV  k 

PGC  =  ^TTTT  C 


Rz(l-k) (l+6k+kz) 


^  Rz (1-k) ^  ^ 


1/2J 


(F-4) 


t 


' 


' 
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Substituting  (F-2)  in  (F-3) 

2/3 


8pkV  4k(o-p)g  . 

ppp  —  _ ^  r _ -I  J- / -3 

R(l+6k+kz)  l3(l-k)yitJ 


( F-5) 


Substituting  in  (3.49) 

..  5k(l-k)2(l-k2)  .4k(o-p)  r,1/3  2/3  1,,.,2  , 

Vb  ‘  4 (l+6k+k2)  [3(T-k7lTj.  1  R  c  +  2(1+k  >  Vc  (F  6> 


For  high  viscosity  liquids,  the  first  term  is  probably  quite  small 
so  that 


Vc  (l+kz)  ^b 


(F-7) 


Substituting  in  (F-5) 


16yV  k 
b 


R(l+kz)  (l+6k+kz)  3  (1-k)  y£V,  ^ 


2k(l+k2) (a-p)  g  1/3 


(F-8) 


PGC 


*5 


Ill 


APPFNPTX  r 


p  mp 77  pponPAM  'PLOT'  PLOTS  ON  Z*Y  SCALF  TPF  COLUMNS 

p  0”  *Z»  MATRIX  AGAINST  TPF  *TPST  COLUMN  D7  TPF 

p  COMMAND  Z  Y  PLOT  X. 


vPT,OT  rrv}u 


m 
r?  i 

r  4  ] 
r  5  ] 
re] 

[7] 

rp] 

m 

[10] 

[11] 

[12] 


V  /  PFrtf  V\C\D\F\G\P\T\J\L\T'%Y\PZ\NP\VT\PT\ST\TSV\V 


\ST\SM 
P£«-' o*°VAn’ 


57V  1  o  f; 


5  10 

■+(  (  0  =  x/(  2 p/1  )  ,pl? )  ,  3  2  1  <ppP) /0  %PL7  tPLl ,PL2 
->PL2  tpP+*)(  2  ,P)p(  tP«-p  ,7?)  ,F 
PP1  :  F«-F[  1  ;  ;  ] 

2  : 7«-l  +  i  (  p  P  )  T  2  ]  - 1 

^<'f[/[/n[;7])-l/!/Pr;7]),r[/Pr;i])-|/F[;i] 

P«- 1  (  2pA  ) *C+ ( C=0 ) *P[1  ;  2  1 ]  +  P[ 1  ;  2  l!  =  0 

F<-(ST[  +  /{  I  0 . 0  001+Pxio*-f7)o  .  5>pr]  )x10*M10»p 

C+(SM*F)x  \  (U/l/P[  ;7]),l  /F[  ;  1  ]  )  xF*SM+l  6  ISM\  1 


4 


[13]  ;i x  n . 5+pr2 ]xpt ;i ]-p[2 ] 

r  1  4  ]  *[;7>!  0.5  +  *[l]xpr;7]-Pri1 

[15]  P+SM* [ ( ( [ / r/P[ ;7] ) , [ /B[ ; 1 ] ) *SM 

[16]  NP+Ol 1  ]  +  ( SMC l]*F[l])x0,iF[l] iSM[ 1 ] 

[17]  7'Z«-G[2]  +  (SW[2]*F[2])x0,iP[2]*S/f[2] 
riPl  0oST+fip~rSv+ix[/+9 

[10]  VL 3  :VT+v/0>NX4-NPxl0*U-STre-fSVl<-T<-l  +  F/  \  10®  |  (NF*0)/N* 

[20]  pt+-\  i  +  io  \pt-i  |p^ip"5  +  (  |/;p)o .  *io*"i+<|>ir/ 

[21]  L+U+ l-(4>(  (C+pW)pl)A.=Pr)i0 

[22]  -*(  (P>7,^-y2,+  [/  7-,  (L+L*T)  ,  (7<0  )x2+F-T)  ,57’[  2  -7SP>S5P [ 
?.-TSV]vL7iU‘VT+r,>T)/  3  2  +126 

[23]  -><  "1  +  126)  ,pSTT  4- PF!']«-r«-l 

[24]  -+PI,  3  ,p/'7P«-£/f[l+~7Fn*“l+i  C 

[25]  P?M -l'7+0[l -7)<t>P7 

[26]  pr«-(  ,P!')xt^,isi(())ppj')pf  ,$r  i*pt  )V .  a(  iP)o  .sii'rrr-i ) ,  (px.p+i- 
r^i/y»+r[  rsn  )pl 

r  2 7  ]  -*(~!'!7’)/2  +  l2R 

[28  ]  P7*r  (  F-  +  /(r,7/)pe7)  +  7/x“l  +  ir]«-llx/7P<0 

[29]  PT+(~( \U+J)e.(  T+J)  ,  i  “l +JW/-7  )  \  ( 1  0  +C,U)pPTtU  pO 

[30]  P7[  iC;  r+7]«-12 

[31]  P7+-’  0123456789“.  »  [1+P2’[  ;  i77-l  ]] 

[32]  ->(~rPF)/  PL  6 

r 3 3 1  7>i  ,/'[2]poxr+pri] 

[34]  P/,4  :/>(  L*PSxC*0  )[1  ,77 [2  ]p0 

[35]  7,[1  +  (P*0)/P[  ;1]>(P*0)//MC  =  P[  ;7])[  .x7 

[36]  -*-(  C* 0  ) /PL  5 

[37]  Wi[0  =  (S/'[2l*2  )  I  0  ,  i/7[2] 

[38]  PLS:PTI  (JpPT)n  ]  #l+r*gW[l])[l  +  0=£.w[l  ]  |P]  ;  ]  ,  (  »  |  »  ,  (pY)pPC 

)  r  !  +  /■;] 


1 

- 

•"  )  l  X 
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a  PROGRAM  'PLOT'  CONTINUED  .... 


*p£onP3<n 

r 3 9 1  •+(  o<r«-o-i  )/P04 

r40]  +  (  U*U+SHZ9l-~TSV+~pfiNP+9P7,)/PL3 

[41]  PL  6  :  ( SM  [2]-9)<M  ,  (  0  0  ,  (  U- 1  )  p  1  ) \PT )  ,  •  ’ 

[42]  +(STl 1  3  2  4  ]  , 1  ) /  1  3  5  7  10  +  12  6 

[43]  'nvrcTl J  APT)  RCALF  rA  CTOP  for  ORDTNATE:  *  ;  5[  1  ]  ,  *Fl  1  ] 

[44]  -+(  0=5?r3]  1/2+126 

[45]  'SCALP  FACTOR  f0p  ORDINATE  \  1  ;1 0*STl  5  ]-l 

[46]  +( 0=57[2 ] )/2+l26 

[47]  *  ORIGIN  AND  SCALP  FACTOR  ”OP  APSCrsSA:  '  ;5[2],*F[ 

2] 

r  4  n  ]  ->f  o=^r[4l  )/o 

[43]  » SCALP  FACTO”  v0p  a  BSC  tSSA  :  '  ;  1  0  *ST[  B  ]  - 1 

[50]  -*>0 

[51]  PLh'TRE  R T  CRT  ARGUMENT  OF  PLOT  MUST  RAVE  RANK  <  3.' 
[51  ]  V 


a  TFE  PROGRAM  ' LTNM'  ESTIMATES  TRE  LEAST  SQUARE  PARAMETERS 
a  ,  VARTANCF  OP  PTTt  RESIDUAL  Errors  AND  CONFIDENCE  LTMITS 
a  n f  PARAMETERS  FOR  A  LTNPAR  MODEL. 


VL  777/f[H]  V 

V  L TNM 1 N ; M ; / ; 595 ; 7 ; C ;SS ; ER ; D 


m 

N+oY 

[2] 

H+pXlli ] 

[3] 

W)+  .  xZ 

[4] 

P+A  t.x(W)  +  .xy 

r  5  ] 

””+-Y-X  +  .  xP 

[B] 

SS+-+  /  PR*  2 

[7] 

V A  R+SS  t  (N  -M) 

[8] 

C+VARxA 

[  9] 

' PARAMETERS  : 

'  \P 

non 

P+(MtM) oO 

rm 

CD  5-<-Mp  n 

[12] 

T+0 

[13] 

LAP:  T+T+ 1 

[  14] 

595[ Tl+2xC[ T\ 7]*0. 

5 

[15] 

0[  r;J>/?[  r  ;7]*“  0 . 5 

riR] 

-*£/!  p  x  \  T <M 

[17] 

C-*-P+  ,  x  C+  .  x  D 

[18] 

'  9  5  P5//7  LIMITS  : 

'  ;  59  5 

[19] 

'VAR  :  '  jF/IP;’  57.0 

DEV  :  ’  ;  IMP*  0.5 

[20] 

’ CORR  MAT  :  •  ;  5 

r  2 1  ] 

'NOR!1  AT,  EPROPS  :  ' 

lERiVAP* 0.5 

r  2  2 1 

WRfl  •  ->•  ^  0 

V 
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p  TFr  PPOCPA!'  '  HOT  •  FSTTMATFS  TPP  I, FAST  SOUAPF  PARAT'FTFPS 

P  , VAPTAllCF  OF  FIT ,  RFSJPUAL  FPPOPS  ART)  CORFTDFNCF  LTPTTS 

r  <9F  PAPAT'.RTFRS  ”0P  A  NON  -L  TNF.AR  MODFL.  TT  FOLLOWS  TPF 

p  TT^/iTvn'r  algorttph  pfvflopfd  py  paro,uappt(rff  .  23). 


V:7f7?rn]V 

7 


•  <7- 


;P2  ;P1  ;7>2 \G\T\TL\  TT \P1  ;P2  ;PP;C 95  ;PPI1  ;PP72  ;7;F7;/ 


r  i  “i  7/1 
r  2 1 

[3]  />py 
r  4  n  tt+ o 

r  5  1  T+(  xff)o  .  =  Xff 

C  G  1  T,T>p.ppT<.v  rjinti  j'H 

r?]  y<-p  psp’<  tp' 

8]  r  ^(^)  +  .x(y-pFr) 

[0]  SS+(K)(  Y-PPT))+ .x(Y-PFT) 

no  1  /«-(  (*?*)+.  XX) 

rill  n-i^np^-T- 

ri2l  r/>o 

[13  1  LAS : TL+TL  +  1 

14D  nir  rL-,TL]<-(A[  7-7  ;  771  +  7/7)*“  0.5 
r  1  5  3  P?.[  r  7  ;7  7  ]•*-(/[  77/ ; 77  1  +  7/. "*//[/)  *”0 . 5 

riRl  +T,AS*\TT,<fl 

f7  1  ni*-(  ni  +  .  x(f?ni  +  .  x  f  /  +  7/”x  n  +  .  xni )  + .  xpi )  +  .  x/7 

ri8l  n2«-(  72+ .  x(f?7)?+  .  x(yl  +  (7/t'7//f/)xr)+  .  x/)2  )  +  .  xP2  )  +  .  *0 

[131  P1+-P+01 

[2  0  1  n  2  *-7  +  72 

r  2 1 1  pp+p 

[22  1  n+-Pi 

[23  1  pn+i^-p  77/ 

[24  1 

[251  PPT2+P  rUCN  LM 

[26  1  51+-(«?(y-PP71  )  )+  .  x  (  Y-PFTl  ) 

[27  1  S?  +  (*)(  Y-PPT?.  )  )+  .  x(  Y-PVP2  ) 

T  2Rl  +?npr*\  (  (  |  (+//  f  7?P-P1  )*F7)  )  )<0 . 0001  ) 

[291  -+0l7Fx\S?.*SS 

3  0(1  -*777  x  i  Pi  >55 

311  [  p<-71 
r  32 1  FS+S1 

[3  3  7  +T,PO 
r  34  1  ORF  :  LAT'+I.A  "*NV 
[35  1  +PNFlx\S2<Sl 
[36]  F+-P1 

637]  SS+S 1 
r  3  a  1 

[39]  />7»P1:P+-P2 

{Ol  5/N-P2 

[41]  -*77(7 

[42]  T!/0:LA!*<-LAP'xRU 

r  4  3  ]  77  7  P 

r  44 1  LnG:P ,ss 

[4  5  1  -*7  PP 

[46]  TPPF : +FOUR* \SS<S1 

[47]  x+f  ps nr  tp' 


' 
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p  ponr-P/ 1M  '  »  COFTTNVFP . 


7wrcn47] 

94  7  ]  7«-P  n«7///f 

4  P  1  r  55+-51 

[49]  rn+.y_pTTri 

fo] 

[51]  POVP:P+-PP 
E2l  FP+Y-PFT 

I!)  3  1  77  r[/^  .  »  p/i  pa  METEPF  ;  »  .  p 

r  s  4  ]  r/  v+ssrp-" 

5B]  n+VAPx  (l?(*?7)  +  .  XX) 

[56]  rns*-" 

[57]  P+T 
r  5  p  ]  r«-n 

[50]  T,ST:T<-T+1 

6D1  r,95rr>2x(C[7;7l*0.5) 

[61]  n[7;  7>5[  P;7]*-0 . 5 

T  6  2 ]  -*>LSTx\T<M 
T  6 3 ]  ' 3 5  PONT  LMT :  ’  ;59  5 

r  64  ]  ’  YA”TANCF  ;  '\VAP\'STD  pry  :  *  ;  STP+VA  P*  0 . 5 

r  6  5  ]  ’5f>Pn  w/2*  :';/)+.xr+.xn 

[66]  ' NO  PM A L  FPPPPS :  ’  \EP*STP 

[67]  ’ MOPFL  L  tKL  tPP  :  '  \L+-VAP  FST  FS 

[67  ]  V 


a  TP  77  PROGRAM  ' FUCN'  STORES  TPF  MODEL  FUNCTIONAL 

o  FORM  FOP  USE  TN  * MOT ' . 


V7i/r//[n]7 

7  PP7VP  P  7/(777  £7? 

rn  ^oupxxip1- 4 

2  1  r  -*TFPEx  \  Lft-3 
3 1  ■+TVf)*\LM=2 

4  I  -+STX*\  LM-6 
5 J  +SEVENx\LM= 7 

r  6  ]  0?1F:PPT+PT1  ]x((  (  (1-7/  )*1  .  5  )  4  77  F )  *  2  x  F  f  2])x((  (7,5-1  )*75)  ) 
r  7  1  ->^7^ 

r  8  1  TJ’O-.PFT+n  n  ]  +  F[2  ]xIVI7**r3  ] 

[9]  +FTVE 

[10]  TFRE'.PFT+P  ri]x(  (75-1  )  *F[  2  ]  )  x  ( 1  -7/1  *2  )  x  ( 1  -XA  *4  ) 

[11]  -*F1TF 

n  21  T*fiiiT>:pvT+p  n  1  x  (  (  (75-1  )  )  *7?r  2  ]x  ( 1  -7/  *2  )  )  x  ( 1  -7/  ) 

[13]  ^rryr 

[14]  5  r7  :  PF  r«-7?  [l]x7’x(*-7?[2-Hl.987xnx(l-*-/?r3]xZ*pr4]) 

[15]  pti/Pj-mO 
V  . 


' 


' 


. 


115 


a  TRF  PROGRAM  '  PSNM'  STORFS  TjrF  DFRTVATTVFS  OF  TRF  MODFL 
a  WTTR  RrSPFOT  TO  TRF  PARAMFTFRS  FOR  USF  7/7  'MOT'. 


vdsnm rnuv 

V  X+P  PSNM  LM i M 
ri]  M+pP 

r  71  X<-(\N)  ol  =  l/f 

r  T  "I  I /■,■*♦=  1 

4i  r  1 7^=7 

5 ir  +TPPFx\LM= 3 

[6]  +FOURx\LM= 4 

[7]  -vSJ/x  iL'*=6 

80  r>wp:^r  ;1  l+P/'rfpri  H 

r  9]  ZT  ;2l<-PPrx2x»(  (  fl-JT/  )*1  .5  )*/7F) 

[10  1  ->7F/7 

[111  7[/r):Z[  ;1  >1 

171  [  tf[  ;2>lMV*ff[3] 

7  1  3  1  Xf  ;3>pr2lx(«F/1K)xI//l/*^r3l 
[14  1  ->7*77 

[15  1  TFRF:X[  ;1  >P/^i7?[ll 

[16]  Z[  ;2>P/77x(»X£-l  ) 

[17]  ->77717 

[18  1  FOVPtXt  ;1  >P/*r*F[l  1 

rim  7r  ;2l>P7,rx(l-^*2)x(#fL9-l  )  ) 

[20  1  ->77/7 

[211  S7T;*[  ;i  >ppr*P[l  ] 

[221  7r  ;2  >-PPrrl  .  987x7 

[23]  ^[;3>p[i]x7x(Z*P[4])x(*-P[2]il.  987x7)x(*-B[3]xZ*P[4]) 

[241  XT  ;4>*[  ;  3 1 xp[ 3 1 x ( «Z ) 

[251  ->7F/7 

[2  61  TpR  :  ->i  0 

V 


' 


